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I .   INTRODUCTION 

Performance  prediction  codes  for  solid  rocket  motors 
model  two  phase  flow  losses  as  functions  of  particle  size. 
In  addition,  particle  size  within  the  grain  port  strongly 
affects  the  damping  of  combustion  pressure  oscillations.   At 
present  these  models  are  based  on  particle  size  data  from 
collected  exhaust  samples  [Ref .  1] .   However,  particle  size 
varies  with  position  in  the  motor  and  other  parameters 
(pressure,  propellant  formulation,  nozzle  design,  etc.) . 
Therefore,  experiments  to  determine  how  particle  size  varies 
in  the  actual  flow  environment  of  the  motor  (i.e.,  across  the 
nozzle)  are  needed  to  validate  the  models  for  two  phase  flow 
losses.   Cramer  [Ref.  2]  and  Karagounis  [Ref.  3]  provide  a 
good  summary  of  the  subject  and  the  Naval  Postgraduate  School 
Combustion  Laboratory  effort  to  obtain  particle  size  data 
across  the  exhaust  nozzle  of  a  solid  rocket  motor. 

The  method  used  in  this  continuing  effort  was  the  dif- 
fractively  scattered  light  technique.   The  diffraction  patterns 
of  light  scattered  by  particles  are  analyzed  to  determine  the 
volume  to  surface  mean  diameter  [Refs.  4  through  11] .   This 
method  has  the  disadvantages  that  size  distributions  cannot 
be  easily  determined  and  particles  larger  than  some  threshold 
size  will  not  be  detected  due  to  the  exceedingly  small  angles 
as  which  they  scatter  light.   However,  it  has  the  advantage 


that  it  is  non-intrusive  ^r.d    in  theory  can  be  used  in  tna 
internal  motor  environment. 

Use  of  this  method  for  particle  sizing  in  solid  rocket 
motors  at  the  Naval  Postgraduate  School  Combustion  Laboratory 
was  begun  by  Karagounis  [Ref .  3] .   The  apparatus  was  sub- 
sequently redesigned  and  the  data  acquisition  equipment  up- 
graded with  the  introduction  of  the  Hewlett  Packard  3054A 
data  acquisition  system  with  an  HP  85  as  the  controlling 
computer  [Ref.  12] .   The  investigation  by  Cramer  and  Hansen 
followed  and  showed  that  propellent  composition  can  limit 
the  application  of  the  technique.   Large  particulate  combustion 
products  in  the  flow  made  particle  size  data  difficult  to 
obtain.   This  was  especially  true  if  only  one  measurement  of 
the  scattering  profile  was  made  during  a  test  firing. 

To  address  this  problem  in  the  present  study  several 
improvements  were  made.   A  cleaner  burning  propellant  was 
obtained  to  reduce  char  agglomerates  in  the  exhaust  products. 
A  more  statistically  valid  data  sample  (multiple  measurements 
during  a  single  test)  was  made  possible  with  added  momory 
in  the  data  acquisition  equipment  and  a  pacing  circuit  which 
allowed  full  use  of  this  memory.   Data  reduction  was  also 
improved  with  a  Hewlett  Packard  98  3  6S  computer  combined  with 
a  more  recently  developed  approach  to  particle  sizing 
presented  by  Buchele  [Ref.  13] .   This  method  is  discussed 
later  in  depth. 
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The  focus  of  this  chesis  "reject  was  the  following; 

(1)  Implement  the  HP  9836S  as  the  system  controller. 

(2)  Expand  the  multiprogrammer  memory  in  order  to  obtain  up 
to  eight  consecutive  scans  of  the  diode  arrays  during 

a  test  firing. 

(3)  Improve  data  reduction  techniques  by  the  method  of 
Buchele  [Ref .  13] . 

(4)  Modify  the  apparatus  and  experimental  procedures  to 
improve  the  angular  resolution  and  to  reduce  extraneous 
light. 

(5)  Certify  the  proper  functioning  and  accuracy  of  the 
apparatus  prior  to  actual  motor  testing. 
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II .   THEORETICAL  BACKGROUND 

A.   GENERAL  DISCUSSION 

The  completely  general  theory  of  scattering  was  developed 
by  Mie  and  is  presented  by  Van  de  Hulst  [Ref .  14] .   The 
light  scattering  characteristics  for  spherical  particles  of 
any  size  are  fully  described  in  a  mathematical  format.   The 
Mie  scattering  functions  contain  Legendre  polynomials  and 
spherical  Bessel  functions  and  fully  treat  the  phenomena  of 
reflection,  refraction,  diffraction,  and  extinction.   The 
full  theory  is  most  often  applied  when  particle  size  is 
approximately  the  same  as  the  wavelength  of  the  incident  light 
Van  de  Hulst  [Ref.  14]  calls  this  the  regime  of  Anomalous 
Scattering . 

For  particle  sizes  much  smaller  than  the  wavelength  of 
light  the  Mie  equations  simplify  to  a  form  which  is  more 
dependent  on  the  index  of  refraction  of  the  particles  and 
less  dependent  on  particle  size.   This  is  called  Rayleigh 
Scattering . 

The  study  of  particle  size  behavior  in  solid  propellant 
rockets  mainly  covers  sizes  much  greater  than  the  wavelength 
of  light.   Scattering  by  large  particles  such  as  these  is 
described  adequately  by  Fraunhofer  diffraction. 
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B.   APPLICATION  TO  LARGER  PARTICLES 

The  ringed  diffraction  pattern  generated  by  a  hole  in 
a  mask,  or  a  number  of  particles  of  the  same  size  is  de- 
scribed by  the  equation: 

2J  (aG)  2 
1(9)  =  [ ± ] 

a9 

where : 

1(6)  describes  the  relative  intensity  of  the  scattered 

light  at  an  angle  theta  (6) 

J,  (a9)  is  the  Bessel  function  of  the  first  kind 

ttD 
a  =  t —   is  the  particle  size  parameter  for  diameter 

D  and  wavelength  of  light  lambda  (A)  . 

Measuring  the  particle  size  for  a  monodispersion  can  be 
accomplished  by  measuring  the  angular  position  of  a  dark  or 
bright  ring  in  the  diffraction  pattern.   For  a  dark  ring  the 
zero  of  the  Bessel  function  corresponding  to  the  ring  is  set 
equal  to  (ct6)  and  particle  size  is  determined  directly  [Ref.  5] 
For  bright  rings  one  sets  (a9)  equal  to  the  corresponding 
maximum  of  the  Bessel  function  and  solves  for  the  diameter. 

The  above  method  is  not  used  for  polydispersions  since 
the  discrete  rings  are  not  observed.   However,  Dobbins,  et  al . 
[Ref.  5]  introduced  a  significant  improvement  in  the  dif- 
fractively  scattered  light  method  of  particle  sizing.   They 
found  that  although  the  method  was  not  directly  able  to 
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determine  distributions  of  sizes,  the  volume  to  surface 
mean  diameter  defined  by 


/°°°N     (D)D3dD 
D         =   -g— L (1) 

JZ  />N     (D)D    dD 

0      r 


where 


N  (D)  is  a  distribution  function  describing  the  pro- 
portion of  particles  with  diameter  (D)  in  the  sample, 
could  be  accurately  measured. 
A  curve  for  sizing  polydispersions  was  presented  which  was 
used  by  Cramer  and  Hansen  [Refs.  2,  12] . 

Two  phase  flow  losses  are  often  calculated  in  terms  of 
D._.   If  the  distribution  of  sizes  in  the  polydispersion  is 
well  behaved  then  D__  and  other  diameters  such  as  D...  can 
be  easily  related  [Ref .  15] .   Reference  5  reported  that  very 
small  particles  in  the  distribution  have  a  minor  influence 
on  the  scattering  profile.   This  makes  the  measurement 
technique  promising  for  the  two  phase  flow  loss  study 
since  very  small  particles  do  not  contribute  significantly 
to  these  losses  and  so  are  of  less  interest. 

Roberts  and  Webb  [Ref.  6]  essentially  confirmed  the 
conclusions  of  Reference  5  and  presented  a  similar  curve  for 
use  in  sizing. 

More  recently,  Buchele  [Ref.  13]  gives  a  good  summary  of 
experimental  techniques  for  particle  sizing  by  measuring 
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diffracted  light.   One  point  of  interest  in  his  report  is 
that  he  represents  the  scattering  profile  of  a  polydispersion 
with  a  function  which  closely  approximates  the  curves  of 
References  5  and  6. 

I  (6)  =  EXP  - ( .57a6) 2 
n 

This  function  from  Reference  13  and  the  curve  from 
Reference  5  were  both  used  in  the  present  study  to  evaluate 
the  apparatus  to  be  used  with  solid  propellant  rocket  motors. 

An  additional  detail  of  measuring  scattering  profiles  is 
covered  by  Van  de  Hulst  [Ref .  14] .   The  wavelength  of  light 
used  in  the  scattering  calculations  depends  on  the  index  of 
refraction  of  the  medium  containing  the  particles.   The 
wavelength  used  in  all  calculations  must  be: 

M 

where : 

An   is  the  wavelength  of  light  in  a  vacuum  and, 
M   is  the  refractive  index  of  the  medium  with  respect 
to  a  vacuum. 
Thus,  the  size  parameter  (a)  becomes: 

ttDM 


a  = 


Ao 


and  the  beam  spread  parameter  becomes: 
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~DM6 

<  r    ^  \        — 

-     7 

0 


Another  consideration  is  presented  by  Gumprecht  and 
Sliepcevich  [Ref .  4] .   Light  scattered  by  particles  in  a 
medium  is  refracted  as  it  crosses  each  interface  of  the 
container  holding  the  medium.   This  is  discussed  further  in 
the  section  on  calibration  and  evaluation  of  the  apparatus. 

Additional  complications  arise  with  the  full  treatment 
of  the  index  of  refraction  of  the  particles  with  respect  to 
the  medium.   But,  for  Fraunhofer  diffraction  alone  this  aspect 
can  be  neglected. 

C.   RESTRICTIONS  AND  SOURCES  OF  ERROR 

Some  restrictions  on  use  of  the  method  are  described  by 
Dobbins,  et  al .  [Ref.  5]  and  were  satisfied  as  described  by 
Cramer  [Ref.  2] .   These  are  related  to  the  size  of  particles, 
the  distance  to  the  detector,  and  some  phenomenon  covered  in 
the  rigorous  Mie  theory. 

One  must  keep  in  mind  also  that  the  curves  developed  for 
polydispersions  are  based  on  the  Upper  Limit  Distribution 
Function  of  Mugele  and  Evans  [Ref.  7] .   This  means  that  no 
particles  with  size  greater  than  approximately  ten  times  the 
mean  should  be  in  the  sample  [Ref.  13] .   This  appears  to  be 
a  mild  restriction.   Van  de  Hulst  [Ref.  14]  describes  the 
criteria  for  single  scattering  and  a  simple  test  to  verify 
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it.   In  general,  as  long  as  the  scattered  intensity  is 
proportional  to  the  number  of  particles  the  mathematics 
remain  simple. 

Sources  of  error  of  the  dif f ractively  scattered  light 
method  are  covered  by  Buchele  [Ref.  13]  and  are  presented 
here . 

(1)  Inaccuracy  of  angular  measurement  or  the  limited  ability 
to  resolve  small  angles  and, 

(2)  Inaccuracy  of  the  intensity  measurement  due  to  extraneous 
light. 

Extraneous  light  includes  all  light  other  than  scattered 
light  from  the  particles-.   Some  examples  are  scattering  from 
an  aperture  or  dirty  test  section  windows.   Refraction  of  the 
beam  due  to  gas  density  gradients  and  image  point  broadening 
from  turbulence  are  others.   Laser  speckle  is  also  extraneous 
light. 

The  sources  of  error  addressed  in  this  investigation  are 
discussed  in  the  related  portions  of  the  paper. 
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III.   EXPERIMENTAL  APPARATUS 

A  photograph  of  the  apparatus  is  presented  in  Figure  1 . 
A  schematic  is  presented  in  Figure  2.   The  light  scattering 
equipment  was  mounted  on  two  optical  benches.   Components  for 
measurements  in  the  exhaust  plane  were  mounted  on  one  bench. 
The  other  bench  held  the  equipment  associated  with  the  motor 
cavity.   The  light  source  was  an  eight  (8)  milliwatt  Helium 
Neon  laser  mounted  on  the  exhaust  bench.   A  collimated  beam 
was  required  so  a  spatial  f ilter/collimator  was  used.   A 
modification  to. this  collimator  is  discussed  later  in  this 
section.   The  collimated  beam  passed  through  a  cube  beam 
splitter  and  the  second  beam  was  diverted  to  a  90  degree 
prism  on  the  other  bench.   The  original  beam  continued  through 
the  motor  exhaust  plane.   The  other  beam  was  routed  through 
the  nitrogen-purged  glass  windows  in  the  motor  housing. 

Each  beam  was  then  intercepted  by  a  physical  stop  located 
in  front  of  its  set  of  receiving  optics.   The  further  the 
stop  was  placed  from  the  test  section,  the  smaller  the  angle 
at  which  scattered  light  could  be  measured.   In  this 
apparatus,  the  stops  were  placed  approximately  30.5 
centimeters  from  the  exit  plane  of  the  test  section.   This 
allowed  a  minimum  angle  of  approximately  .008  radians  to  be 
measured.   Light  scattered  at  angles  greater  than  this  was 
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not  intercepted  and  continued  past  the  edges  cf  the  stop. 
The  stop  served  to  keep  the  transmitted  beam  out  of  the 
measuring  optics  and  thus  reduce  extraneous  light.   The 
stops  also  improved  optics  alignment.   This  is  discussed 
under  calibration  and  evaluation  of  the  apparatus. 

The  scattered  light  passed  through  a  narrow  pass  filter 
which  admitted  only  light  of  the  Helium  Neon  frequency. 
This  filter  served  to  reduce  extraneous  light  from  the  external 
surroundings . 

An  objective  lens  of  50  centimeter  focal  length  was 
located  behind  the  narrow  pass  filter.   This  lens  imaged 
onto  a  photodiode  array  the  scattering  profile  of  the  particles 
in  the  test  section.   The  shadow  of  the  beam  stop  was  also 
imaged  since  the  stop  was  between  the  test  section  and  the 
objective  lens.   This  was  a  limitation  which  is  discussed 
under  calibration  and  evaluation. 

The  photodiode  arrays  were  the  same  units  used  by  Cramer 
and  Hansen  [Refs.  2,12].   Each  array  contained  1024  silicon 
photodiodes  on  a  single  chip  with  25  micron  spacing.   The 
accompanying  circuits  provided  a  sampled  and  held  output 
which  was  essentially  analog  except  for  switching  transients. 
At  the  end  of  each  diode  scan  there  was  a  delay  before  the  next 
scan.   During  this  delay  the  diodes  were  reset  and  allowed 
to  measure  the  intensity  of  the  scattering  profile  again. 
The  scanning  of  the  diode  array  repeated  continuously.   The 
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actual  sampling  time  of  the  array  was  about  3  4  milliseconds 
with  a  delay  between  scans  of  about  6  milliseconds. 

The  50  centimeter  focul  length  of  the  objective  lens 
combined  with  the  dimensions  of  the  diode  array  provided  a 
half  angle  field  of  view  of  about  3  degrees  for  mediums  of 
refractive  index  near  one.   The  effective  field  of  view  was 
reduced  to  about  2.3  degrees  for  calibrations  when  the 
refractive  index  of  a  Plexiglas  container  and  water  was 
taken  into  account . 

The  laser  beam  collimator  mentioned  previously  at 
first  produced  a  beam  one  centimeter  in  diameter.   A  lens 
in  the  collimator  was  changed  to  reduce  the  beam  diameter  for 
several  reasons.   Extraneous  light  would  be  generated  if 
a  large  beam  impinged  on  the  aperture  of  the  motor  test 
section  window.   Also,  if  the  aft  beam  was  larger  than  the 
motor  exhaust  jet  it  would  be  refracted  in  the  density 
(and  refractive  index)  gradient  between  the  exhaust  and  air. 

The  last  part  of  the  apparatus  was  the  rocket  motor 
itself.   It  was  the  one  used  by  Cramer  and  Hansen  [Refs.  2,12] 
and  in  the  present  study  served  only  for  aligning  the  optics. 
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IV.   DATA  ACQUISITION  SYSTEM 

A.   NEW  CONTROLLER 

Hansen  [Ref.  12]  describes  the  major  components  of  the 
Hewlett  Packard  3054A  data  acquisition  system.   A  list  of  the 
manuals  relevant  to  this  study  is  in  Table  1.   The  HP  85 
computer  used  by  Hansen  and  Cramer  was  replaced  with  an 
HP  9836S  as  system  controller.   This  newer  computer  has  far 
more  capability  than  the  HP  85,  including  a  choice  of  more 
powerful  operating  systems.   The  system  used  for  this 
study  was  Basic  Extended  2.1. 

The  data  acquisition  program  written  by  Hansen  needed 
minor  modification  to  acquire  multiple  consecutive  scans  of 
the  photodiodes.   Some  different  I/O  commands  such  as  those 
which  transfer  data  to  the  disk  were  also  incorporated. 
The  revised  version  of  this  program  is  listed  in  Appendix  A. 
A  general  flow  chart  is  presented  in  Figure  3 . 

The  9836S  has  two  internal  disk  drives  which  were  used 
to  store  the  data  after  acquisition.   The  data  from  both 
diode  arrays  was  stored  in  the  same  file.   The  eight  (8) 
scans  of  the  motor  cavity  were  first,  followed  by  the  four  (4) 
exhaust  scans. 
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B.  MODIFICATIONS 

The  memory  capacity  of  the  Mul tiprogrammer  unit  was 
increased  so  that  multiple  consecutive  scans  of  the  diode 
arrays  could  be  recorded  during  a  motor  firing.   This  would 
provide  a  more  statistically  valid  measurement  of  particle 
size.   Fluctuations  of  scattered  light  intensity  for  a 
polydispersion  need  to  be  integrated  over  time  or  averaged 
to  provide  a  more  appropriate  measurement. 

In  order  to  fully  use  the  memory  added  and  make  data 
management  easier  the  data  acquisition  system  needed  to  be 
modified.   The  memory  consisted  of  three  (3)  cards,  each 
with  a  capacity  of  4096  values.   The  fact  that  this  was  a 
multiple  of  1024  (the  exact  number  of  photodiodes)  meant 
that  the  idle  period  between  scans  needed  to  be  excluded 
from  the  data.   If  this  was  not  done,  one  (1)  less  scan  per 
card  would  have  been  acquired  and  locating  the  scans  in 
the  overall  block  of  memory  would  have  been  more  difficult. 

It  was  also  necessary  to  chain  two  of  the  memory  cards 
together  in  a  way  which  would  allow  one  card  to  be  filled 
and  then  the  other.   A  schematic  of  the  data  acquisition 
system  is  presented  in  Figure  4. 

C.  PACING  AND  MEMORY  CONTROL  CIRCUITS 

The  timing  clock  and  blanking  pulse  of  the  photodiode 
circuitry  provided  the  means  for  pacing  data  acquisition. 
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Specific  results  of  the  modifications  were: 

(1)  Memory  space  was  fully  utilized  and  management  of  the 
multiple  scans  made  easier. 

(2)  A/D  conversions  of  the  data  were  made  exactly  when  a 
diode's  output  was  on  line  and  steady.   Thus,  the  analog 
filter  used  in  the  previous  study  was  no  longer 

needed  to  suppress  the  switching  spikes  on  the  data  line. 

The  following  is  a  description  of  the  signals  and  circuits 
used  to  modify  the  data  acquisition  system.   All  voltage 
levels  were  TTL .   A  timing  diagram  in  Figure  5  shows  the 
relations  between  signals.   A  schematic  of  the  circuit  is 
presented  in  Figure  6 . 

The  clock  pulse  was  a  positive  going  spike  at  a  frequency 
of  about  30  KHz.   This  clock  controlled  all  circuits  of 
the  photodiodes.   It  ran  continuously,  even  during  the  blank 
period  between  scans  when  the  diode  output  was  clamped  at 
zero  volts. 

The  blanking  pulse  was  a  signal  which  fell  to  zero  at 
the  beginning  of  each  scan.  It  then  went  positive  at  the 
end  of  the  scan  and  remained  high  until  the  next  scan  began. 

The  clock  pulse  was  used  to  drive  a  pulse  shaper  (mono- 
stable  multivibrator) .   This  ensured  that  the  voltage  levels 
through  the  rest  of  the  circuit  would  not  accidentally  fall 
below  the  TTL  threshold.   The  pulse  width  of  the  shaper 
was  adjusted  so  that  the  negative  going  edge  of  each  pulse 
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would  occur  after  the  switching  transient  on  the  data  line 
had  decayed.   This  negative  going  edge  would  eventually 
trigger  the  A/D  converter  to  store  the  output  of  each  diode. 

The  blanking  pulse  was  inverted  and  connected  to  an  AND 
gate  along  with  the  pulse  shaper  output.   The  output  of  this 
AND  gate  is  shown  in  Figure  5  as  the  pulse  shaper  signal  held 
low  between  scans  of  the  diode  array.   This  was  the  basic 
signal  which  paced  data  acquisition. 

This  basic  trigger  signal  was  connected  to  an  AND  gate 
along  with  the  output  of  the  Multiprogrammers  Timer  Pacer 
card.   In  this  way  the  trigger  would  not  reach  the  A/D  until 
the  Timer  Pacer  output  a  pulse.   This  enabling  pulse  from 
the  Timer  Pacer  was  at  least  as  long  as  the  time  for  eight 
scans  of  a  diode  array. 

The  controller  programed  the  Timer  Pacer  to  produce  the 
pulse  when  the  Timer  Pacer  received  a  trigger  from  the  blanking 
pulse.   In  this  way,  data  acquisition  began  at  the  start  of 
a  scan  and  no  data  was  taken  during  the  time  between  scans. 

The  circuit  to  chain  the  memory  cards  together  was 
basically  an  OR  gate  used  as  a  negative  logic  AND  gate.   The 


end  of  conversion  (EOC)  signal  of  the  A/D  and  the  (FULL)  signal 
of  one  memory  card  were  connected  to  the  gate.   When  both 
signals  went  low  the  second  card  was  then  able  to  store  data 
from  the  A/D.   This  arrangement  is  shown  in  the  schematic  of 
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the  data  acquisition  system  in  Figure  4.   The  automatic  lock- 
out feature  of  the  memory  cards  when  full,  and  the  relatively 
slow  rate  of  data  throughput  made  it  unnecessary  to  control 
other  handshake  lines  [memory  card  manual,  Table  1] . 

The  circuit  was  designed  to  handle  four  (4)  memory  cards, 
so  no  modification  will  be  necessary  if  one  more  card  is 
added  to  the  system.   This  would  provide  an  additional  four 
scans  of  the  exhaust  beam. 
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V.   DATA  REDUCTION 

The  data  reduction  programs  written  by  Hansen  were  not 
used  for  this  study.   The  new  computer  lended  itself  to 
another  approach.   Its  memory  capacity  made  it  unnecessary  to 
chain  programs  together  and  polynomial  curve  fitting  was 
eliminated  in  favor  of  interactive  graphics.   Avoiding 
polynomial  fits  preserved  the  nature  of  the  raw  data  so  that 
one  had  a  better  feel  for  the  parameters.   The  data  reduction 
program  "RDC"  is  listed  in  Appendix  A.   Figure  7  is  a  general 
flow  chart  for  the  program. 

The  following  is  a  description  of  the  program.   The  user 
was  first  prompted  for  values  needed  to  analyze  a  given 
data  set.   For  example,  the  wavelength  of  the  laser  used  and 
the  index  of  refraction  of  the  medium  must  be  known  for 
any  data  set.   Next,  one  had  the  choice  of  reducing  raw 
data  scans  or  reviewing  a  reduced  data  file.   For  raw 
data  one  chose  either  the  exhaust  or  motor  cavity 
beam  data. 

Raw  data  was  plotted  on  the  CRT  and  any  obviously  erroneous 
scans  were  excluded  from  further  reduction.   The  valid  scans 
were  averaged  to  obtain  a  mean  scattering  profile.   The  mean 
intensity  profile  taken  before  particles  were  introduced 
was  then  subtracted  from  that  taken.,  with  particles  present. 
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This  corrected  for  the  characteristics  of  individual  photo- 
diodes  and  extraneous  light  which  was  independent  of  the 
particles . 

A  symmetric  moving-average-type  of  digital  filter  was 
then  applied  to  the  profile  to  achieve  some  smoothing. 
This  type  of  digital  filter  was  chosen  for  simplicity  and 
because  it  does  not  have  the  phase  lag  of  analog  filters 
[Ref .  16] .   Preserving  the  phase  of  the  data  was  necessary 
to  retain  angular  resolution.   Another  advantage  of  filtering 
in  the  software  rather  than  hardware  was  that  raw  data 
files  remained  unmolested. 

The  scattering  profile  was  then  analyzed  using  interactive 
graphics.   If  earlier,  one  chosen  to  review  a  reduced  file, 
program  execution  began  here. 

One  had  to  normalize  a  scattering  profile  in  order  to 
compare  it  to  the  theoretical  curves  for  polydispersions . 
The  scattered  intensity  on  the  centerline  of  the  beam  was 
the  correct  value  to  use  for  normalization  but  was  unmeasurable 
due  to  the  beam's  presence. 

The  other  unknown  was,  of  course,  the  particle  size.   These 
two  variables  (centerline  intensity  for  the  measured  profile 
and  D-.„  for  the  theoretical  profile  of  normalized  intensity 
vs.  (9,))  were  adjusted  using  interactive  graphics  until 
the  curve  for  polydispersions  coincided  with  the  data.   In 
this  way  the  mean  diameter  of  particles  was  determined. 
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The  second  reduction  technique  used  was  the  direct 
application  of  the  method  presented  by  Buchele  [Ref .  13]  . 
The  equation  for  the  polydispersion  curve: 

I  (9)  =  EXP  -  (  .57  a6) 2 

n 

was  applied  at  two  points  of  the  scattering  profile.   This 
gave : 

12/11  =  EXP  -D2[(82  -  82)(.57  -^/X)2] 

Solving  this  for  the  diameter  gave: 


2    2     2 
,2/(62 


d  =  [-l  (i2/ii)  (x/.57TT)V(e;  -  e;-)]': 


The  computer  would  sweep  through  the  data  using  many 
values  of  thetal  (9,)  along  with  several  angle  ratios  to 
determine  theta2  (8~) .   The  results  were  presented  graphically 
as  particle  size  vs.  thetal  for  each  angle  ratio  (8  /9  )  . 

In  actual  practice  the  range  of  useable  angles  depends  on 
the  appratus,  and  the  quality  of  the  data.   Therefore,  in  order 
to  interpret  the  results  one  must  have  previously  inspected 
the  data.   The  interactive  graphics  routine  was  well  suited 
to  this  and  provided  a  hard  copy  for  inspection. 

After  reducing  a  set  of  raw  data  the  mean  scattering 
profile  was  stored  on  disk  for  later  review. 
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VI.   CALIBRATION  AND  EVALUATION 

A.   IMPROVEMENTS 

The  geometry  of  the  apparatus  used  in  the  investigation 
by  Cramer  and  Hansen  is  compared  with  that  of  this  study 
in  Figure  8.   In  the  previous  study  the  transmitted  beam 
was  allowed  to  enter  the  receiving  optics.   The  beam  was 
focused  off  the  diode  array  a  few  millimeters  from  the  first 
diode.   This  was  necessary  to  avoid  damaging  the  diodes 
but  introduced  some  uncertainty  in  angle  measurements. 
The  intense  image  of  the  beam  along  with  scattered  light  from 
the  receiving  optics  produced  a  high  level  of  extraneous 
light.   In  the  present  study,  stops  were  used  to  intercept 
the  beam  before  reaching  the  receiving  optics.   These  stops 
provided  several  advantages.   A  high  intensity  beam  could 
be  used  while  producing  little  extraneous  light.   Optics 
alignment  was  also  improved.   This  reduced  error  in  angle 
measurement.   Alignment  was  accomplished  using  a  neutral 
density  filter  to  reduce  beam  intensity  and  protect  the 
diodes.   A  schematic  of  the  apparatus  is  in  Figure  2.   The 
laser,  collimator,  beam  splitter  and  prism  were  positioned  so 
that  the  beams  passed  through  the  appropriate  measurement 
areas.   The  narrow  pass  filter  and  imaging  lens  were  then 
positioned  so  that  the  beam  entered  on  the  centerline.   The 
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photodiode  array  was  then  moved  using  a  three-axis  micrometer 
so  that  the  focused  beam  fell  on  the  first  diode.   The  beam 
stops  were  then  put  in  place  and  the  neutral  density  filter 
removed.   In  this  way  measurements  commenced  exactly  from 
the  optical  axis  of  the  beam. 

Procedures  were  also  refined  to  account  for  the  bending 
of  light  rays  as  they  passed  through  the  walls  of  the  particle 
container.   As  noted  earlier,  the  index  of  refraction  of  the 
container  and  the  medium  containing  particles  affects 
scattering  measurements.   A  Plexiglas  box  held  the  particle 
samples  and  a  magnetic  stirrer  kept  the  samples  suspended 
in  water.   The  index  of  refraction  of  the  Plexiglas  and 
water  combination  was  measured  using  a  simple  technique. 
A  microscope  was  used  to  measure  the  ratio  of  actual  depth 
to  apparent  depth  for  Plexiglas  and  water.   The  index  of 
refraction  was  determined  to  be  1.39.   This  value  was 
applied  to  the  data  to  convert  the  measured  scattering  profile 
to  that  actually  produced  in  the  medium  containing  the 
particles . 

B.   RESULTS 

Calibration  results  are  summarized  in  Table  2.   Initial 
tests  were  done  with  two  samples  wihch  were  basically  mono- 
dispersions  of  large  particles.   Figure  9  shows  the  measured 
profile  of  scattered  light  for  glass  spheres  ranging  from 
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37  to  44  microns  in  diameter.   This  profile  was  obtained  by 
placing  the  focussed  beam  just  far  enough  from  the  first 
diode  to  avoid  saturation  with  no  particles  present.   The 
diodes  located  at  angles  less  than  about  .01  radians 
saturated.   The  first  bright  ring  for  particles  of  about  40 
micron  diameter  was  visible  near  .02  radians.   Figure  10 
shows  a  profile  for  the  same  particles,  illustrating  use  of 
the  beam  stops  to  avoid  diode  saturation  and  improve  angle 
measurements.   In  this  case,  the  first  diode  was  located 
exactly  on  the  centerline  of  the  beam  as  discussed  above. 

Results  for  a  sample  of  53  to  63  micron  glass  spheres  are 
shown  in  Figure  11.   The  center  lobe  was  nearly  completely 
missed  but  the  first  two  bright  rings  were  seen  near  .014 
and  .022  radians.   The  first  two  dark  rings  near  .01  and 
.019  radians  were  also  seen.   The  method  described  earlier 
of  setting  the  beam  spread  parameter  (a6)  equal  to  the 
zeroes  and  maximums  of  J, (a8)  was  used  to  calculate  a  size 
of  about  58  microns.   Also  shown  in  Figure  11  is  the 
theoretical  profile  for  a  polydispersion  with  D_„  =  54 
microns . 

Various  polydispersions  of  either  glass  spheres  or 
aluminum  oxide  powder  were  then  tested.   These  polydispersions 
consisted  of  fairly  large  particles.   Results  are  shown  in 
Figures  12  through  17.   These  tests  showed  that  the  apparatus 
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had  two  distinct  modes  of  operation.   If  the  particle  concen- 
tration was  very  high,  or  if  large  particles  dominated  the 
polydispersion,  many  of  the  diodes  at  the  smaller  angles 
would  saturate.   This  left  only  the  data  at  larger  angles 
useable.   When  many  diodes  saturated,  the  theoretical  curve 
given  by  Dobbins,  et  al.,  was  used  to  determine  size.   This 
was  done  because  this  curve  was  valid  for  the  larger  angles 
and  lower  relative  intensities.   The  curve  from  Buchele 
[Ref.  13]  was  not  valid  for  values  of  the  beam  spread 
parameter  greater  than  three  (3)  . 

For  low  particle  concentrations  and/or  small  particles  the 
data  proved  more  accurate  at  the  smaller  angles.   If  no 
diodes  were  seen  to  saturate  then  one  knew  the  measurement 
was  in  the  higher  intensity  part  of  the  center  lobe.   Here 
the  curve  given  by  Buchele  was  quite  satisfactory  for  sizing. 

The  smallest  particles  tested  were  five,  ten,  and  twenty 
micron  polystyrene  spheres.   The  bright  rings  for  these 
particles  occurred  at  angles  too  large  for  the  apparatus  to 
measure.   For  these  samples  the  diodes  did  not  saturate. 
Both  the  Gaussian  curve  fit  and  the  two  angle  method  were 
used  to  obtain  D__.   These  results  were  especially  consistent. 
It  should  be  noted  that  the  two-angle  method  uses  the  equation 
for  the  Gaussian.   If  the  measured  profile  matched  the 
Gaussian  exactly,  then  D  _  would  be  the  same  for  any  (6,) 
and  angle  ratio   (9„/9-)  employed.   Some  variations  in 
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calculated  D    due  to  the  imperfect  fit  are  obvious  in 
Figures  19,  21  and  23. 

A  scanning  electron  microscope  was  used  to  photograph 
the  types  of  particles  tested.   These  photographs  are  shown 
in  Figures  24  through  27.   Equation  (1)  was  used  along  with 
these  photographs  to  calculate  some  of  the  values  of  D  ~ 
in  Table  II.   Calculations  of  D  9  for  the  polystyrene  were 
arrived  at  using  the  manufacturers  data  on  size  distributions 
The  photographs  generally  confirm  the  validity  of  the 
technique . 
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VII .   CONCLUSIONS  AND  RECOMMENDATIONS 

The  results  of  the  calibration  tests  showed  that  the 
apparatus  is  capable  of  accurately  measuring  mean  particle 
size  for  a  broad  range  of  mean  diameters.   It  was  found  that 
the  technique  was  most  accurate  if  the  theoretical  profile 
fit  or  the  two-angle  method  were  applied  at  the  smallest 
possible  scattering  angles. 

The  rocket  exhaust  is  likely  to  attenuate  the  beam 
somewhat,  reducing  the  problems  related  to  diode  saturation 
at  small  scattering  angles.   Thus,  measurements  should  be 
possible  using  the  high  intensity  part  of  the  center  lobe. 
This  should  make  data  reduction  less  ambiguous.   Actual 
testing  should  begin  with  measurements  at  the  exhaust 
plane  of  the  motor.   These  should  be  compared  with  collected 
exhaust  samples  to  validate  the  use  of  the  apparatus  in 
an  actual  motor  environment.   Measurements  in  the  motor 
cavity  would  then  be  interpreted  based  on  the  correlation 
between  exhaust  samples  and  exhaust  measurements. 

It  is  also  recommended  that  the  index  of  refraction  of 
the  combustion  gases  be  investigated.   A  literature  search 
for  an  estimate  of  the  index  of  refraction  would  probably 
be  satisfactory. 
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TABLE  I 
ELECTRONICS  MANUALS 

1.  HP  Memory  Cards  Model  6970B  Operating  Manual 

2.  HP  Timer  Pacer  Card  Model  69737A  Operating  Manual 

3.  HP  Analog  to  Digital  Converter  Card  Model  69736A 
Operating  Manual 

4.  HP  Users  Guide,  "Using  the  9826  and  9836  Computers  with 
the  6942A  Multiprogrammer" 

5.  Basic  Language  Reference  Guide  with  Extensions  2.0  for 
Series  200  Computers 


35 


CO 

E-i 

-J 

r> 

00 

w 
cc 

z 

o 

I — I 
< 

i— l 
< 


w 

cq 
< 
E- 


00 

«fr 

LO 

h    c 

o     o 

LO 

[•O 

vo 

O 

O 

H        H 

• 

• 

• 

• 

4-> 

4-J 

H            U 

cn 

CV1 

tH      -H 

O 

fO 

e 

4-> 

SO  c 

C    <U 

CO 

o 

■Hg            en 

LO 

tn 

HO)            C 

(U     H       N     O 

LO 

cr-i 

o 

o 

+->   3     ro  ?h 

. 

• 

4-> 

+-> 

+->   oo  Q      u 

■=3- 

C-^ 

r— 1 

o 

CO 

LO 

oS   o3           -h 

C"J 

^3" 

•=* 

CO 

CNI 

^J- 

U   0)            E 

LO 

CN 

oo  S 



rH  X) 

v— '  CD 

* 

•K 

■K 

+-)               00 

* 

■K 

•K 

•K 

•K 

•K 

<r 

<r 

C    03             C 

CN 

CNI 

O    r— 1        N     O 

t*- 

rvi 

\0 

0) 

0) 

•  H    3      Kl    (h 

• 

• 

• 

• 

. 

• 

0)       • 

CJ 

■H    U  Q        U 

** 

o 

rH 

CO 

■st 

LO 

oo    ao 

en    oo 

03   H               -H 

rH 

oo 

fO 

LO 

OJ 

•H 

•H 

3   o3           E 

&H 

En 

cru 

w 

CD 

O 

CD 

H" 

ro 

r» 

rH              00 

vO 

H 

tn 

«* 

O 

CO 

O    CD    C 

LO 

O 

•HMO 

O 

o 

o 

o 

O 

o 

CN 

LO 

•H   -H     !_ 

H 

4-> 

■M 

+-) 

4-> 

4-J 

H  OO    u 

21 

21 

03           -H 

hO 

>o 

LO 

r~~ 

ro 

rH 

Oh            E 

rH 

ro 

LO 

cu 

CD 

CU 

0)        rH 

e 

C 

£ 

E 

E 

i— I       03 

0) 

CD 

CD 

3 

3 

U       -H 

H 

H 

H 

CO 

oo 

oo 

C    0) 

C     0) 

•H          $H 

X 

>^ 

^ 

oo 

oo 

oo 

•H  Td 

•H   X3 

•H        CD 

4-> 

+J 

+J 

03 

03 

03 

E  -h 

E  -h 

H        4-> 

1/5 

!/) 

oo 

H 

rH 

rH 

3   X 

3    X 

03       03 

X 

>x 

X 

C3 

CD 

03 

rH    O 

rH     O 

Ph        S 

i— 1 
O 
Pn 

rH 
O 

rH 
O 
PL, 

< 

< 

00 

CD 

H 


U 


O 

H 


■K 


oo 
O 

4-> 

o 
a. 

00 

g 
o 

H 


36 


Figure  1.   Photographs  of  Light  Scattering  Apparatus 
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APPENDIX  A 

PROGRAM  LISTINGS 

10  i**********        ACQDTA  :  FOR  9336S  ****************** 

20  !**********  ACQUIRES  MULTIPLE  SCANS  CF  BOTH  DIODE  ARRAYS  ***************** 

30  !**********       AND  STORES  THEM  ON  DISK         ***************** 

40  !**********  BERT  HANSEN,  KEN  GSAHAM, KELLY  HARRIS  1984   ***************** 

50  OPTION  BASE  1 

60  COM  Ql*[123,O2$[12],Tl$[20],Dl$[12],D2S[12],Address(4),E(4096) 

70  COM  A(1024)  BUFFER, B(4096)  BUFFER, C(4096)  BUFFER  ,D(40?6)  BUFER 

80  OUTPUT  709;  "AR"  'ANALOG  RESET 

90  ASSIGN  gHultilO  TO  72310  ICLEARS  THE  UAXE-UP  SERVICE  REQUEST 

100  ENTER  §Multil0;Gql,Sq2,8q3,Qq4,Qq5)G;q6         !CF  THE  HULTIPROGRAMMER 

110  MAT  Address=  (3) 

120  ENABLE  INTR  7 

130  CLEAR  722 

140  PRINT  USING  '3/" 

150  PRINT  "ENTER  THE  FILENAMES  OF  THE  DATA  FILES  TO  BE  CREATED  (e.q.  RAUlfRAW2 

)" 

160  PRINT  ■  12  CHARACTERS  MAXIMUM,  EACH' 

170  PRINT  USING  7' 

180  PRINT  ■         AN  EMPTY  DISK  MUST  BE  IN  THE  LEFT  DISK  DRIVE  ■ 

190  PRINT  '* 

200  PRINT  '  THE  DATA  FILES  UILL  NEARLY  FILL  A  DISK' 

210  PRINT  USING  '/" 

220  MASS  STORAGE  IS  ': INTERNAL, 4,0"        !  CHANGE  THIS  LATER  IF  NECESSARY 

230  INPUT  '         INPUT  FILE  NAMES  NOW  -  <FILENAME1,FILEXAME2>  ",D1$,D2$ 

240  M5=1024 

250  Zz$=": INTERNALS,  1"  !  STRING  INDICATES  MASS  STORAGE 

260  CREATE  BDAT  Dl$4Zz$,6144,16   !  12  SCANS  CF  1024  g  2*8  BYTES  PER  RECORD 

270  CREATE  BDAT  D2${.Zz$,6144,16   !  1024*12/2=NUHBER  OF  REC0RDS=6144) 

280  PRINT  USING  "?" 

298  PRINT  "      DATA  UILL  BE  STORED  ON  DISKETTE  WITH  FOLLOWING  FILE  NAMES:' 

300  PRINT  USING  '///' 

318  PRINT  "  NO  PARTICLES  FILENAME  =  ";DU 

320  PRINT  B  PARTICLES FILENAME  =  ";D2$ 

338  PRINT  USING  '//' 

340  PRINT  '  IS  THIS  A  CALIBRATION  ?   ENTER  '  Y  '  IF  YES  ' 

330  PRINT  u  ANYTHING  ELSE  IF  NO  ' 

360  INPUT  R$ 

378  PRINT  USING  '//' 

380  PRINT  '  BE  SURE  LASER  IS  ON' 

390  PRINT  "  PRESS  [CONTINUE]  WHEN  READY" 

400  PAUSE 
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410 
423 
430 
440 
450 
460 
470 
480 
490 
500 
510 
520 
530 
540 
553 
560 
570 
580 
590 
600 
610 
620 
638 
640 
650 
660 
670 
680 
690 
708 
718 
720 
738 
748 
750 
760 
770 
780 
790 
808 
810 
820 
838 
840 
850 
868 
870 
880 
890 
908 


TC 


R$= 


«V«  T 


nt' 


H  470 


PRINT  'USING  "2" 

INPUT  "ENTER  THE  THRESHOLD  PRESSURE  TO  TRIGGER  THE  DEVICES  (psi)',S8 

INPUT  "ENTER  TIME  DELAY  FROM  THRESHOLD  PRESSURE  (s8c)\T8 

!        TIMES  ARE  INTERPRETED  BY  THE  COMPUTER  IN  SECONDS  DOWN  TO  .001 

V0=S8/151.5  !  CALIBRATION  OF  PRESSURE  TRANSDUCER 

D3i=Dl$  !  NO-PARTICLES  STRING  NAME 

GOSUB  Multiprog 

GOSUB  Storedata 

IF  R$="Y"  THEN  1800 

!************************************************************************ 

Contcheck:  (Continuity  checking 

OUTPUT  7C9;*AC20"  "CONNECT  3456A  DVM  TO  IGNITION  FIRING  LINE 

LOCAL  LOCKOUT  7  (KEEPS  OPR.  FROM  SWITCHING  FRONT  PANEL  TO  LOCAL 

R7=225  !  IN  THOUSANDS  OF  OHMS 


OUTPUT  722; "HSM032SW2SO1L1S0F4R1T3SX1 " 


'RESISTANCE  MEASUREMENT 


H  RESETS  DVM;SM002  SETS  SERVICE  REQUEST  MASK  WHERE  002  IS  OCTAL  REP,  OF 
THE  SERIAL  POLL  MASK  BYTE;SW2  TELLS  WHICH  TERMINAL  SWITCH  IS  USED;S01 
SYSTEM  OUTPUT  MODE  ON  -  WAITS  FOR  CONTROLLER  TO  HANDSHAKE ;L1  LOAD 
INTERNAL  MEMORY  ON;S0-FUNCTION  SHIFT  CFF;F4-TWO  WIRE  CONNECTION  TO  DVM 
R1-AUT0RANGING;T3-SINGLE  TRIGGER jQ-LOAD  INTERNAL  MEMORY  OFF;Xl -EXECUTE 
STORED  PROGRAM. 

GOSUB  Reading 

RS=V 

IF  R8<R7  THEN  GOTO  Ok 
!************<*********CONTINUITY  CHECK********************************* 

PRINT  USING  '?" 

PRINT  "CONTINUITY   CHECK   BAD!!!!' 

PRINT  "RECHECK  BEFORE  PROCEEDING,  WHEN  CHECKED,  PRESS  [CONTINUE]" 

BEEP  3000,. 3 

WAIT  .1 

BEEP  100,1,0 

PAUSE 

GOTO  Contcheck 
OkiPRINT  USING  '8' 

DISP  ■ 

PRINT  ' 

OUTPUT  709;"AC21" 

PRINT  ' 

PRINT  USING  "/■ 

PRINT  "BE  SURE  VISICORDER  IS  SET  UP  TO  RUN  ON  PROPER  SCALE  WITH  LAMP  ON." 

PRINT  USING  "/" 

OUTPUT  722; "HSM002SW2Z0SO1L1FL0 , 01STIS0F1R1T3Q 

DISP  "  STANDING  BY  FOR  IGNITION" 

PRINT  "  STANDING  BY  FOR  IGNITION" 

BEEP  2000,. 1 

OUTPUT  722} "XI ' 

GOSUB  Reading 

R9=ABS(V) 

IF  R9< 10  THEN  GOTO  870 


CONTINUITY      CHECK      IS       O.K." 
BE  SURE  NITROGEN  IS  ON* 
! CONNECT  DVM  FOR  IGNITOR  VOLTAGE  MEASUREMENT 
DVM  CONNECTED  TO  FIRE  SWITCH" 


VOLTAGE  ON  IGNITER 


'VOLTMETER  TRIGGER 
•READS  VOLTMETER 


!  12  VOLTS  ON  IGNITOR 
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910 

920 

930 

940 

950 

968 

970 

980 

998 

1000 

1010 

1023 

1030 

1040 

1050 

1060 

1070 

1080 

1090 

1180 

1110 

1120 

1130 

1145 

1150 

1160 

1170 

1180 

1190 

1280 

1210 

1220 

1230 

1240 


EHPUTER  CLGCX  TO  ]ET  ELAPSED  TI-.E3 
i  CONNECT  PRESSURE  XDUCER  TO  DYM 
!  TRIGGER  VOLTMETER 


THRESHOLD  PRESSURE 
!  TIME  DELAY 


O0=TIME£ATE  ;  USES 

OUTPUT  709;"AC22' 

OUTPUT  722; 'XI' 

GOSUB  Reading 

R9=ABS(V) 

IF  R9<V0  THEN  GOTO  930 

WAIT  T8 

01=TIHEDATE 

IF  R${)"Y"  THEN  1060  ! IF  ACTUAL  RUN  THEN  SKIP  SOME  LINES 

!  FOR  CALIBRATIONS  INTRODUCE  PARTICLES  AND  THEN  CONTINUE 

PRINT  CHR$(12) 

PRINT  USING  '////////' 

PRINT  '  INTRODUCE  PARTICLES  THEN  PRESS  KEY  *  9  TO  TAKE  DATA" 

ON  KEY  9  LABEL  '  TAKE  DATA*  GOTO  1060 
Standby:  GOTO  Standby 

GOSUB  Hultiprog  !  TAKE  PARTICLE  DATA 

LOCAL  7 

D3$=D2i 

GOSUB  Storedata 

PRINT  " 

PRINT  ' 

PRINT  ' 

GOTO  End 
Hultiprog;  ! 

OUTPUT  723 

OUTPUT  723 

OUTPUT  723 

OUTPUT  723 

OUTPUT  723 

OUTPUT  723 

OUTPUT  723 

OUTPUT  723 

OUTPUT  723 


! PARTICLE  DATA  FILE  NAME 


ELAPSED  TIME: 
ELAPSED  TIME: 


FIRE  TO  MULTIPROGRAHMER  CALL  = 
FIRE  TO  MEMORY  INTERRUPT 


-,01-00 
'-.02-00 


DATA  STORED  ON  DISK  WITH  FILENAMES  (\D1$;*)  AND  (")D2$;")* 


25C 


1 
1260 


"00,2,3,12,13^ 

■CC, 5, 6,9,101'  ICLEARS  THE  ARM, BUSY  AND  EOP  OF  MEMORY  CARDS 
"CC,1,4,11T"   .'CLEARS  SAME  FOR  A  TO  D'S  AND  TIMER  PACER 
■SF, 2, 3,1, .001,12,1"     !  THE  (1)  IS  2'S  COMPLIMENT  BINARY 
^,5,3,1,  .001, 12, T"   !  THE  .001  IS  THE  LEAST  SIGNIFICANT  BIT 
'SF,9,3,1,.001,12,T'   !  THE  12  IS  FOR  12  BIT  WORD  SIZE 
■SF,12>3>1 , .001 >12,T"  !  SINCE  THE  A  TO  D  IS  12  BIT 
■yF,3,1023T,WF,6,4095T' 
"UF,10 ,4095T,WF, 13,40?5T" 
!  SETS  REFERENCE  WORD  FOR  WHEN  TO  STOP  TAKING  DATA  AND  GENERATE  INTERRUPT 
OUTPUT  723 j  *WF , 2 . 1 , 1 T , WF , 5 . 1 , 1 T , WF , 9 . 1 , 1 T , UF ,  1 2 . 1 , 1 T "     » SETS  FIFO  MODE 
OUTPUT  723;'WF,3,l,0T,UF,3.2,OT,WF,3.3,0T,WF,13.1,OT,UF,13.2,0T,yF,13.3,OT 


1270  OUTPUT  723;'WF,6.1,OT,WF,6.2,0T,yF,6.3,OT,WF,10.1,0T,WF,10.2,0T,yF,10.3)0T 


1280         !  SETS  COUNTER  AND  POINTERS  OF  2ND  MEMORY  CARD  IN  EACH  PAIR  TO  0 

1290  OUTPUT  723;,AC,3T,AC,6T,AC,10T,AC)13T,!ARMS  CARDS  WHICH  GENERATE  INTERRUPT 

1300  OUTPUT  723;-UF,4.2>QT,WF,4>lST"  !  TIMER  PACER  GIVES  1  PULSE  OF  1  SEC  WHEN 

1310         !  TRIGGERED  BY  THE  BLANKING  PULSE  (PLENTY  OF  TIME  FOR  8  SCANS) 

1328  WAIT  2 

1330  K=SP0LL(723)  IUAIT  FOR  HEM  INTERRUPT 

1340  IF  K064  THEN  GOTO  1330 

1350  OUTPUT  723;'yF,4.2,lT"    !  MAY  BE  UNNECESSARY  TO  ALTER  TIMER  PACER  SINCE 

1360  !         MEMORY  CARDS  HAVE  AUTOMATIC  LOCKOUT  BUT  FOR  NOU  WE  WILL  DO  IT 

1370  02=TIMEDATE  " 

1380  ON  ERROR  GOTO  Err  trap       !  NEEDED  TO  READ  ARMED  CARD  INTERRUPT  LIST 

1390  SEND  7;UNL  MLA  TALK  23  SEC  12   !  SPECIFICALLY  ASKS  FOR  INTERRUPT  LIST 
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1400 
1410 
1420 
1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 


ENTER 

PRINT 


7;Addressl*) 


var_read 
Mencards 

PRINT  " 

MAT  Address=  (0) 

OFF  ERROR 

OUTPUT  723;,DC,3,6,13(10T' 

OUTPUT  723; "HI ,2 , 1 024T " 

ENTER  72305;A(») 

OUTPUT  723;'HI,5,4096T* 

ENTER  72305;B(») 

PRINT  " 

OUTPUT  723;,HI,12,4096T1 

ENTER  72305;C(«) 

OUTPUT  723;'MI,9,4096T* 

ENTER  723C5;D(*) 

ENABLE  INTR  7;8 

PRINT  * 

(-1)*A 

(-1XB 

(-1)*C 

(-1>*D 


READ  WHICH  CARDS  INTERRUPTED 


MEMORY  CARDS  WHICH  GENERATED  INTERRUPTS 
SLOTS*  =  ';Address(*> 


tfE 


(DISSARH  MEM  CARDS 

'  SET  UP  CARD  TO  BE  READ 

GETS  DATA  FROM  1024  MEMORY  BOARD 


EXHAUST 


!  GETS  DATA  ROM  4096  MEMORY 
DATA  ENTERED' 


BOARD 


MAT  A= 
MAT  B= 
HAT  C= 
HAT  D= 

RETURN 


MOTOR    DATA  ENTERED' 
THE  1024  CARD  IS  INCLUDED  BUT  NOT  SAVED.  IT  DIDN'T 
PERFORM  WELL.  COULD  BE  REPLACED  BY  A  4096  CARD. 


DIODE  VOLTAGES  ARE  NEGATIVE  SO  SIGNS  ARE  CHANGED 


1620  Storedata;! 

1630  ASSIGN  gDiskfile  TO  D3$iZz$ 

1640  ASSIGN  gfiuffl  TO  BUFFER  C(») 

1650  ASSIGN  §Buff2  TO  BUFFER  D(») 

1660  ASSIGN  ?Buff3  TO  BUFFER  B(«) 

1670  CONTROL  8Buffl ,3; 1 ,32763,1 

1680  CONTROL  8Buff2, 3; 1,32767,1 

1690  CONTROL  8B«ff3, 3)1,32767,1 

1700  TRANSFER  BBoffl  TO  SDiskfile 

1710  WAIT  FOR  EOT  Miskfile 

1720  TRANSFER  0Buff2  TO  gDiskfile 

1730  WAIT  FOR  EOT  BDiskfile 

1740  TRANSFER  gBuff3  TO  EDiskfile 

1750  UAIT  FOR  EOT  SDiskfile 

1760  ASSIGN  gDiskfile  TO  * 

1770  assign  mm  TO  * 

1780  ASSIGN  gBuff2  TO  * 

1790  ASSIGN  §Boff3  TO  * 
1800 


!  MOTOR  CAVITY  4  SCANS 

!  MOTOR  CAVITY  4  SCANS 

!  EXHAUST     4  SCANS 

.'SETS  BUFER  POINTERS  TO  FULL 

(INTERFACE  REGISTERS  SECTION  OF 

(LANGUAGE  MANUAL 

!  ORDER  OF  DATA  ON  THE  DISK  IS 

!  MOTOR  CAVITY— 8  SCANS 

!  EXHAUST      -  4  SCANS 

(WAIT  BECAUSE  OVERLAPPING 

(TRANSFERS  ARE  NOT  WANTED 

(CLOSE  I  /  O  PATHS 


ASSIGN  BBuffl  TO 

ASSIGN  gBuff2  TO 

ASSIGN  gBoff3  TO 

RETURN 
1810  Reading:! 
1820  STATUS  7,1 }AB 
1830  ENTER  722;V 
1840  ENABLE  INTR  7;8 
1850  RETURN 

1860  Err_trap:  IF  ERRN=159  AND  ERRL(Var  read)  THEN  Hencards 
1870  PRINT  ERRM$  (EVEN  IF  THE  ERROR  WAS  NOT  THE  ONE  PLANNED 

1880  -GOTO  hencards  (FOR  PROGRAM  EXECUTION  CONTINUES 

1890  End:  END 


(CHECKING  STATUS  BEFORE  READING 
(VOLTMETER  IS  A  FORMALITY  TO 
(CLEAR  THE  SERVICE  REQUEST 
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10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

111 

120 

130 

140 

150 

160 

R 

170 

180 

198 

200 

210 

228 

230 

ON' 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 


i ******************************************************************* 


******************* 
******************* 
******************* 
******************* 


RDC 
PLOTS  RAW  DATA 

AVERAGES 
FILTERS 


*********  *********** 
******************** 
******************** 
******************** 
******************** 
******************** 


*******************  DETERMINES  MEAN  DIAMETER 

*******************  BY  INTERACTIVE  GRAPHICS 

*******************    AND  THE  TWO-ANGLE  METHOD    ******************** 

*******************   Robert  Kellv  Harris    ******************** 

*******************         1984         ******************** 

******************************************************************* 

OPTION  BASE  1 

COM  /Hrdaaus/  Av2(1024) 

COM  /Gauss/  T1(1024),G(1024) ,L 

COM  /Max/  M7,M5,Xt,Yt,XM,Yn 

COM  /Readata/  B, P,H,Q3$[201, 9451201 ,Zz$l2fl],YH8192)  BUFFER ,Y2( 81 92)  BUFFE 


COM  /Two/  Avl(1024),M,Ml,F 
DIM  Scans(8),X(1024) 

COM  /Plots/  P1$[203,P2$[20],P3$[20],P4$[20] 
INTEGER  Graf (1:12480)  BUFFET? 
Choose:  PRINT  CHR$(12) 
Old=0 
PRINT  ■  TO  LOOK  AT  ANY  DATA  FILE  THE  PROGRAM  NEEDS  SOME  STARTING  INFORMATI 


PRINT  " 

PRINT  ■ 

PRINT  " 

PRINT  " 

PRINT  "• 

PRINT  ' 

PRINT  " 

PRINT  " 

BEEP 

INPUT  ' 

PRINT  CHR$(12) 

PRINT  "   TO  ACCOUNT  FOR  REFRACTION  OF  LIGHT  AT  BOUNDARIES  BETWEEN' 

PRINT  " 

PRINT  '  THE  MEDIUM  AND  AIR  YOU  MUST  ENTER  THE  INDEX  OF  REFRACTION 

PRINT  "" 


TO  ACCOUNT  FOR  THE  CHANGE  IN  WAVELENGTH  IN  THE  MEDIUM' 

ENTER  THE  INDEX  OF  REFRACTION  OF  THE  MEDIUM' 

WATER  =  1.33" 
AIR  =  1.8" 
ESTIMATE  OF  EXHAUST  =1,1" 

THIS  VALUE  ADJUSTS  THE  COMPARISON  CURVE  TO  THE  MEDIUM'. Ml 
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3^0 
400 
410 
420 
430 
440 
450 
460 
470 
480 
490 
500 
510 
520 
530 
540 
550 
560 
570 
580 
590 
600 
618 
620 
630 
640 
650 
660 
670 
680 
690 
700 
710 
720 
730 
740 
750 
768 
778 
788 
790 
800 
810 
828 
830 
840 
850 
860 


GF  THE  COMBINATION  OF  THE  MEDIUM  AND  ITS  BOUNDARY1 
THIS  VALUE  IS  APPLIED  DIRECTLY  TO  THE  DATA" 

AIR=1.0  PLEXIGLASS  &  WATER  =  1.39" 


PRINT  ' 

PRINT  " 

PRINT  ' 

PRINT  " 

PRINT  ' 

PRINT  " 

PRINT  "ESTIMATE  OF  MOTOR  CAVITY  COMBUSTION  PRODUCTS  &  WINDOW  =  1.22" 

PRINT  ,9 

PRINT  '  ESTIMATE  OF  ROCKET  EXHAUST  =  1.1  OR  1,0' 

INPUT  M 

PRINT  CHR$(12) 

INPUT  '  ENTER  LASER  WAVELENGTH  (HeNe=. 6328 ,Ar=, 488) ',L 

INPUT  'ENTER  THE  DESIRED  PLOTTING  INTERVAL  OF  DIODE  ARRAY  (2>4,6)',H 


!ttn  FOCUL  LENGTH  OF  OBJECTIVE  LENS 
!ah  DIODE  SPACING  MAY  BE  .03 


F=460 
Diod=.825 

FOR  1=0  TO  1023 

X(I+1)=I+1 

TKI+l)  =  (I*Diod)/(F«H) 

NEXT  I 
FRINT  CHR$(12) 
CALL  Displavl(Old) 
IF  Old  THEN  CALL  Review(AvlU)) 
IF  Old  THEN  GOTO  Gauss 
BEP 

INPUT  ■  INPUT  TWO  FILENAMES  NOW  (NO-PART  ,PART)")S3$1Q4$ 

INPUT  ■      ENTER  '!'  FOR  MOTOR     '2'  FOR  EXHAUST* ,P 


CREATE  AN  ARRAY  OF  DIODE  NUMBERS 
CREATE  AN  ARRAY  OF  THETA 


!JUST  PRINTS  INITIAL  REMARKS  ON  CRT 
ITHESE  TWO  LINES  ALLOW  FOR 
! REVIEWING  REDUCED  DATA  FILES 


Zi$=':INTERNAL,4,r 
CALL  Readata 
CALL  Shift 
BEEP 
Screen:  ! 
PRINT  CHR$(12) 
PU='VOLTAGE* 
P2*='DI0DE  NUMBER" 
P3i=  "VOLTAGE  vs.  DIODE' 
P4$='RAW  DATA' 
M5=1024 
Xt=64 
YT=.04 
Xn=4 
Y«=5 

CALL  Plot(Y2(*)) 
GSTORE  Craft*) 
CALL  Dataplot(B,YK*),H) 
PRINT  ' 
PRINT  " 
PRINT  * 
INPUT  J 


!  USES  LEFT  DISK  DRIVE  FOR  RAW  DATA  FILES 


ICLEARS  SCREEN 

!  STRINGS  FOR  PLOTS 

!  SET  UP  VALUES  FOR  PLOTS 


!  DRAWS  AXES  ,ETC... 
!  STORES  GRAPHICS  DISLAY  JUST  MADE 
!PLOT  NO-PARTICLES  DATA 
NO-PARTICLES  PLOTTED" 

HOW  MANY  SCANS  SEEM  TO  BE  GOOD  ?' 
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373  PRINT  " 

385  PRINT  '        WHICH  SCANS  ARE  GOOD  ?..  .ie. . .1 ,2,4, 5,7, INCLUDE  LAST  CO 
H.MA"  !LAST  COMMA  IS  REQUIRED  OR  YCU  HAVE  TO  HIT  CONT. .  TWICE 

890  INPUT  ScansO) 

900  PRINT  CHRSU2) 

910  PRINT  USING  '/////////' 

920  PRINT  '  AVERAGING  THE  SELECTED  SCANS1 

930  GRAPHICS  OFF 

940  CALL  Averaqe(J,Scans(*>,Yl(*),AvlU)) 

950  PRINT  USING  'g' 

968  GLOAD  Graft*)     ILOADS  GRAPHICS  ARRAY  RATHER  THAN  WASTE  TIME  RE-DRAWING 

970  GRAPHICS  ON 

988  CALL  Result(Avl(*),X(*),H)  !  PLOTS  1024  ELEMENT  ARRAYS 

990  PRINT  USING  '//' 

1000  PRINT  '  Average  Intensity  No-Particles' 

1010  ON  KEY  0  LABEL  "AVERAGE1"  GOTO  Screen 

1020  ON  KEY  1  LABEL  'PLOT-PARTICLES'  GOTO  1060 

1030  PRINT  '  PRESS  KEY  *  0  TO  RE-AVERAGE  OR  *  1  TO  CONTINUE" 

1040  Standby:  !  MANUAL  CALLS  THIS  INTERRUPT  DRIVEN  PROGRAMMING 

1050     GOTO  Standby  !  LOOPS,  WAITING  FOR  USER  TO  DECIDE 

1060  PRINT  USING  "?' 

1070  OFF  KEY  1  IHELPS  AVOID  CONFUSION  BY  CLEARING  THAT  BOX 

1080  CALL  Dataplot(B,Y2(*),H)  ! PLOT  PARTICLES  DATA 

1090  PRINT  USING  7" 

1100  PRINT  '  PARTICLE  DATA  PLOTTED' 

1110  PRINT  '       FOR  A  HARD  COPY  OF  THIS  RAW  DATA  PRESS  KEY  *  6" 

1120  PRINT  '  ■ 

1130  PRINT  ■  OR   TO  CONTINUE  PRESS  KEY  *  1" 

1140  ON  KEY  6  LABEL  "HARD  &  RAW"  GOTO  Raw 

1150  ON  KEY  1  LABEL  'CONTINUE  "  GOTO  Select 

1160  GOTO  Standby 

1170  Raw:  CALL  Plot(Y2U)  ,1)         !  THE  ONE  (1)  IS  AN  OPTIONAL  PARAMETER 
1180     CALL  Dataplot(B,Y2(*),H)      !  WHICH  IS  USED  TO  GET  HARD  COPIES 

1190  Select:! 
1200     OFF  KEY  6 
1210     OFF  KEY  1 

1220  PRINT  CHR$(12) 

1230  PRINT  USING  '////' 

1240  PRINT  '  HOW  MANY  SCANS  SEEM  TO  BE  GOOD  ?" 

1250  INPUT  J 

1260  PRINT  '  WHICH  SCANS  ARE  GOOD  ?...ie. ..1,2, 4,5, 7, INCLUDE 

LAST  COMMA" 

1270  INPUT  ScansU) 

1280  PRINT  CHR$(12) 

1290  PRINT  USING  "/////////' 

1300  GRAPHICS  OFF 

1310  PRINT  "  AVERAGING  THE  SELECTED  SCANS" 

1320  CALL  Averaged, Scans(*),Y2(*),Av2(*)) 

1330  PRINT  USING  "g* 

1340  P4t="AVERACED  SCANS' 
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1350  CALL  Plot(Av2(*)) 

1360  GRAPHICS  ON 

1370  GSTORE  Graft*) 

1380  CALL  Resalt(Av2(«),X(*),H) 

1390  PRINT  USING  '//' 

1400  PRINT  ' 

1410  ON  KEY  2  LAEEL  "AUERAGE2"  GOTO  1080 
ON  KEY  3  LABEL  "SUBTRACT"  GOTO  1450 


Averaae  Scattered  Intensity' 


1420 

1430  PRINT  " 
CONTINUE" 

1440  GOTO  Standby 

1450  PRINT  USING  "g" 

1460  OFF  KEY  3 

1470  HAT  Avl=  Av2-Avl 

1480  GLOAD  GrafU) 

1490  CALL  Result(Avl(*),X(*),H) 

1500  PRINT  USING  "///" 

1510  PRINT  " 
o-Particles" 

1520  PRINT  " 

1530  PRINT  " 

1540  ON  KEY  6  LABEL 

1550  ON  KEY  1  LABEL 

1560  GOTO  Standby 

1570  CALL  PlotCAvl (*), 1 > 

1580  CALL  Result(A«l(*),X(*))H) 

1590  PLOTTER  IS  3, "INTERNAL" 

1600  PRINT  USING  "//" 

1610  PRINT  " 

1620  PRINT  ' 

1630  PRINT  " 

1640  PRINT  " 

1650  INPUT  " 

1660  IF  Fil=0  THEN  Gauss 

1670  P4$="FILTERED  DATA" 

1680  CALL  Filter(Avl(*),Fil) 

1690  PRINT  CHR$(12) 

1700  CALL  Pbt(Avl<*)) 

1718  IF  Old  THEN  HAT  X=  Tl 

1720  CALL  ResBlt(Avl(*),X(* 

1730  PRINT  ' 
o-Particles" 

1740  PRINT  ' 
TER" 
1750 
1760 
1770 
1780 
1790 
1800 
1810 


HIT  KEY  *  2  TO  RE-AVERAGE  OR  KEY  t  3  TO 


(SUBTRACTS  NO-PARTICLES  FROM  PARTICLES 


Plot  of  the  Difference  Between  Particles  and  N 


KEY  *  6 
KEY  *  1 


FOR  HARD  COPY" 
TO  CONTINUE  " 


HARD  AVERAGE"  GOTO  1570 
FILTER  "  GOTO  1600 


ENTER  THE  NUMBER  OF  TIMES  YOU  WISH  TO  APPLY" 
THE  DIGITAL  FILTER  FOR  SMOOTHING" 
EXAMPLE********  10" 
TAKES  ABOUT  1.5  MINUTES" 
ZERO  IS  ALSO  ACCEPTABLE", Fil 


H) 


PRINT  " 
PRINT  ' 
PRINT  ' 
PRINT  " 
PRINT  " 

ON  KEY  8  LABEL 
ON  KEY  9  LABEL 


Plot  of  the  Difference  Between  Particles  and  N 

AFTER  APPLICATION  OF  A  DIGITAL  FIL 

for  a  HARD  COPY  * 
PREPARE  the  PLOTTER  and" 
PRESS  KEY  *  8" 
OR" 
PRESS  KEY  *  9  TO  NORMALIZE' 


■HARD  FILTERED"  GOTO  1838 
'  NORMALIZE  '  GOTO  Gauss 
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1820  GOTO  Standby 

1830  CALL  Plot(Avl(*).n                    !  The  1  is  for  hard  copy 

1840  CALL  Res»lt(Avl(*),X(*),H) 

1850  OFF  KEY  8 

I860  Gauss:  OFF  KEY  2                    'KEY  1  OPTION  IS  STILL  VALID 

1870  OFF  KEY  8 

1880  ON  KEY  1  LABEL  '   FILTER  '  GOTO  1600  !  ALLOWS  ONE  TO  SMOOTH  OLD  DATA 

1890  !  WHEN  REVIEWING  REDUCED  FILES 

1900  PRINT  CHR$(12) 

1910  !Av2  Array  is  Normalized  BUT  Avl(»)  is  Still  Saved  For  Re-work  if  Needed 

1929  MAT  Av2=  Awl   !  TO  BEGIN  WITH  THE  ARRAY  IS  ASSUMED  TO  BE  NORMALIZED 

1930  !  AND  THE  USER  CHANGES  THIS  WITH  INTERACTIVE  GRAPHICS 
1940  !  IN  THE  SUBROUTINE  "CoMoare* 

1950  Pl$="NORMALIZED  INTENSITY" 

1960  P2J='THETA  (rad)"                          !STRINGS  FOR  PLOTS 

1970  P3$="INTENSITY  vs.THETA" 

1980  P4$='CURVE  FIT  RESULTS" 

1990  M5=.05                               !SET  UP  UALUFS 

2000  Xt=,004                                IFOR  PLOTS 

2010  Yt=.l 

2920  X«=2 

2030  Yn=2 

2040  CALL  Plot(Av2(*)) 

2050  CALL  Result(Av2(*),TK*),H)      IPLOT  OF  NORMALIZED  INTENSITY  vs,  THETA 

2060  GSTORE  GrafU)               !  SAVES  THE  SCREEN  IMAGE 

2070  CALL  Co«pare(Avl(*),H,H, Ml, Graft*)) 

2080  GRAPHICS  OFF 

2090  ON  KEY  4  LABEL  'OTHER  ARRAY"  GOTO  Choose 

2100  ON  KEY  5  LABEL  "  TWO-ANGLE"  GOTO  Buchele 

2110  ON  KEY  9  LABEL  '  NORMALIZE"  GOTO  Gauss 

2120  ON  KEY  6  LABEL  "STORE  DATA"  GOTO  Storedata 

2130  ON  KEY  7  LABEL  '   QUIT   *  GOTO  Quit 

2140  CALL  Menuel                                                      !  PRINTS  OPTIONS  ON  THE  SCREEN 

2150  GOTO  Standby 

2168  Buchele:     ! 

2170  CALL  Twoanqle 

2180  CALL  Menuel 

2190  GOTO  Standby 

2208  Storedata:! 

2218  CALL  Store 

2220  CALL  Menuel 

2230  GOTO  Standby 

2248  Quit:  ! 

2250  END 

2260  SUB  Averaged, Scans(*),Y(*),Av(*))  (AVERAGES  SELECTED  SCANS 

2270  MAT  Av=  (0)         UNITIALIZES  THE  ARRAY  LOCAL  TO  THIS  ROUTINE 

2280  FOR  1=1  TO  J  STEP  1 


2290        K=(Scans(I)-l)«1024+l 
2300        L=Scans(I)*1024 
2310  FOR  I1=K  TO  L 


THIS  COUNTER  IS  THE  BEGINNING 
AND  THIS  ONE  THE  END  OF  BLOCKS 
OF  1024  IN  THE  OVERALL  DATA 
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2320 
2330 
2340 
2350 
2360 
2378 
2380 
2390 
2480 
2418 
2420 
2430 
2448 
2450 
2460 
2470 
2480 
2490 
2500 
2510 
2520 
2530 
2540 
2550 
2540 
2570 
2580 
2590 
2600 
2610 
2620 
2638 
2648 
2650 
2660 
2670 
2680 
2690 
2780 
2718 
2720 
2730 
2740 
2750 
2760 
2770 
2780 
2790 
2800 
281 


12=1 1-X+l 

Av(I2)=Av(I2)+Y(Il) 

NEXT  II 
NEXT  I 
MAT  Av=  Av/tJ) 
SUBEND 
SUB  Dataplot(B,Y(*),H) 
COH  /Max/  H7,H5,Xt,Yt,Xfl,Yn 
LDIR  0 
LORG  4 

FDR  J=l  TO  B 

MOVE  500+J*50,M7-,05 

LABEL  J; 

LINE  TYPE  1 

K=(M)*1024+1 

L=J*1024 

MOVE  1,Y(K) 

FOR  I=K  TO  L  STEP  H 
X=I-K+1 
PLOT  X,Y(I) 
NEXT  I 
NEXT  J 
PENUP 
SUBEND 

SUB  Result(Iav(*),D(*),H) 

CLIP  ON 

MOVE  D(l),Iav(l) 

FOR  1=1  TO  1024  STEP  H 
DRAW  D(I),Iag(I) 
NEXT  I 
PENUP 
SUBEND 
SUB  Plot(Y(»), OPTIONAL  Hard) 
COM  /Max/  H7;M5,XT,Yt,X«,Y« 


!  THIS  COUNTER  13  ALWAYS  BETWEEN 
!  1  AND  1024 


i*****THE  AVERAGE*****! 

IB  IS  4  OR  8(THE  NUMBER  OF  SCANS) 


'EACH  SCAN 
!  SEE  NEXT  LINE 
HELPS  KEEP  TRACK  OF  EACH  SCAN  AS  IT  APPEARS 

IBEGINNING  OF  EACH  SCAN  AND 
!THE  END  WITHIN  THE  TOTAL  BLOCK 
!  MOVE  TO  THE  FIRST  POINT 

!  THIS  GIVES  1  TO  1924  FOR  ABSCISSA 


!THE  AVERAGE  INTENSITY  IS  PLOTTED 


!IF  THE  OPTIONAL(Hard)IS  RECEIVED 
!THE  FIGURE  GOES  TO  THE  PLOTTER 


COM  /Plots/  P1$[20],P2*[20],P3*[20],P4$[20] 


M7=MAX(Y(*))*1.1 
M6=,1*M7 
PRINT  " 
GINIT 

SELECT  NPAR 
CASE  1 

PLOTTER  IS  3, "INTERNAL" 
CASE  2 

PLOTTER  IS  705,'HPGL" 
END  SELECT 
GRAPHICS  ON 
CSI2E  5, .6 


!  A  SCALING  VARIABLE 
!  ANOTHER  SCALING  VARIABLE 

IGINIT  IS  JUST  GOOD  PRACTICE  SO  YOU  KNOW 

•WHERE  YOU  ARE  BEGINNING 
'THIS  DETECTS  IF  THE  HARD  COPY  IS  DESIRED 
!IT  COULD  BE  DONE  WITH  IF  THEN  LOGIC 
!BUT  IS  PRESENTED  FOR  FAMILIARIZATION 
•SINCE  IT  IS  MORE  POWERFUL  IN  COMPLEX 
(SITUATIONS 

!  TO  BE  ABLE  TO  SEE  THE  PLOT 


0  LDIR  0 
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ocnn  i  0w~  * 

2930  MOVE  75,95 

2840  LABEL  USING  "K";P4S 

2853  HOVE  75,90 

2860  LABEL  P3$ 

2870  MOVE  75,20 

2880  LABEL  P2$ 

2895  LDIR  90 

2900  HOVE  18,60 

2910  LABEL  PI? 

2920  VIEWPORT  30,125,30,85  !  SCREEN  UNITS  FOR  MARGINS 

2938  FRAME  !  DRAWS  A  BOX 

2940  WINDOW  0,H5,-H6,M7  !  EXTENT  OF  X  AND  Y 

2950  AXES  Xt,H6,0,0,X*,l,5 

2960  LDIR  9 

2970  CSIZE  3,. 5 

2980  LORG  8 

2998  CLIP  OFF 

30B0     FOR  I=-H6  TO  H7  STEP  H7/10  (NUMBER  THE  Y  AXIS 

3010    HOVE  0,1 

3020     LABEL  USING  ,*,HD.DDDD';I 

3030     NEXT  1 

3040  CSIZE  3,. 6 

3050  LDIR  90 

3060    LORG  8 

3070     FOR  1=0  TO  H5  STEP  Xt  INUHEER  THE  X  AXIS 

3088     HOVE  I,-H6 

3090    LABEL  USING  ,*,K"}I 

3109     NEXT  I 

3119  PENUP 

3120  SUBEND 

3130  SUB  Cortpare(Avl(*),H,H, HI, INTEGER  Graf(*)) 

3140  COH  /Gauss/  T1(*),G(*),L 

3150  COH  /Hrdqaus/  Av2<») 

3160  COH  /Hax/  H7,H5,Xt,Yt,Xrt,.Y« 

3170  BIH  T(l:1024),At(l:17),Dobbins(l:17),Tld(l:17) 

3180  DATA  3.,3.5,4.,4,5,5.,5.5,6.,6.5,7.,7,5,8.,8,5,9.,?,5,18.,10,5,11. 

3190  DATA  .0556,. 035,. 0206,. 014,. 0106,, 0081,. 00605,, 005,. 00374,. 0032,, 00248, 

.0022, .00185, .0016, .00135, .0012, .001 

3200  READ  At(») 

3210  READ  Dobbins(») 

3228       OFF  KEY  'GETS  RID  OF  ALL  LABELS  ON  KEYS 

3230       ON  KEY  3  LABEL  ,HENUE1  GOTO  Subexit 

3240       PLOTTER  IS  3,'INTERNAL"     !  IN  CASE  A  HARD  COPY  WAS  JUST  HADE 

3250       Hard=0  !  (0)  SO  ONE  DOESN'T  EXIT  TOO  SOON 

3260       Centerlme=l  !  THE  INITIAL  NORMALIZING  VALUE 

3270       D=20.  !  INITIAL  PATICLE  MEAN  DIAHETER  IN  MICRONS 

3280  Change:  CSIZE  4,. 6 

3290       PRINT  USING  *//////' 


75 


3300 

3310 

3320 

3330 

3340 

3350 

3360 

3370 

3380 

3390 

3400 

VE 

3410 

3420 

3430 

3440 

3450 

3460 

3470 

SIZE" 

3480 

3490 

3500 

3510 

3528 

3530 

3540 

3550 

3560 

3578 


PRINT  ' 

MOVE  0,1 

CLIP  ON 

Con=PI*D*Ml/L 

FOR  1=1  TO  1024  STEP  H 

T(I)=TKI)*Con 

G(I)=EXP(-((.57*T(I))A2): 

IF  T(I))3.0  THEN  3400 

DRAW  T1(I),G(I) 

NEXT  I 

MAT  Tld=  At/tCon) 


FOR  1=1  TO  17 

DRAW  TldCI).Dobbins(I) 

NEXT  I 

PENUP 
IF  Hard  THEN  SUBEXIT 
PRINT  USING  ■/////////* 

PRINT  ' 


STAND  BY  FOR  CURVE" 


!A  CONSTANT 

!Theta-bar  for  sone  Dbartnote  Ml) 
IGaussian  curve  for  this  Dbar 

IGAUSSIAN  NO  GOOD  FOR  Theta-bar)3 
IGaussian  carve  vs,    theta 

ITHETA  FOR  POINTS  ON  DOBBINS'CUR 


.'EXITS  ROUTINE  IF  A  HARD  COPY  WAS  JUST  MADE 


USE  THE  KNOB  TO  VARY  THE  PARTICLE 


PRINT  ' 

PRINT  '  ' 

PRINT  " 

PRINT  " 

PRINT  '  " 

PRINT  ' 

PRINT  USING  '//////' 

ON  KEY  6  LABEL  'HARD  COPY"  GOTO  Hardgauss 

ON  KEY  9  LABEL  "NORMALIZE"  GOTO  Divide 

ON  KNOB  .05  GOTO  Pulse 

GOTO  Wait 

PRINT  USING  "8" 

Strngl$=*Mean  Diameter  " 

Strna2$="  Microns' 
Count=KNOBX 
D=DR0UND(D+Ceunt/15>2) 
GLOAD  Graft*) 
LDIR  8 
LOSG  8 

MOVE  M5,.8*H7 

LABEL  StrngmVALt(DHStrnq2* 
GOTO  Change 


3580  Wait; 

3590  Pulse 

3600 

3610 

3628 

3630 

3640 

3650 

3660 

3670 

3680 

3690 

3700  Hardaauss:! 

3718    PRINT  ' 

3720    PRINT  " 

3730    PRINT  * 

3740    PAUSE 

3750    Hard=l 

3760    CALL  Plot(Av2(*),l) 

3770    CALL  Result(Av2(»)  ;TK*),H) 


OR  HIT  KEY  *  6  FOR  HARD  COPY" 

KEY  *  9  ALLOWS  YOU  TO' 
RE-NORMALIZE" 

HIT  KEY  *  3  TO  GET  OUT' 


!(.05)  IS  TIME  INTERVAL  IN  WHICH 

!  PULSES  FROM  THE  KNOB  ARE 

!  COUNTED  AND  THIS  NUMBER  IS 

!  USED  BY  THE  INTER R ACTIVE 

!  GRAPHICS  TO  VARY  THE  PARTICLE 

!  SIZE  AND  PLOT  THE  ASSOCIATED 

!  GAUSSIAN  APPROXIMATIONS  OF 

!  SCATTERING  PROFILES 


PREPARE  the  PLOTTER" 
PRESS  CONTINUE  fir" 
HARD  COPY" 

!  SO  THAT  SUBEXIT  OCCURS  AFTER  HARD  COPY 
!  1«*HARD  COPY 
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3780  PRINT  USING  "?" 

3790  LORG  8 

3808  MOVE  H5,.8«H7 

3810  LDIR  0 

3820  LABEL  Strnql$4VAL«(D)4Strng2S 

3830  MOVE  0,1 

3840  GOTO  Chanqe 

3850  Divide:! 

3860  PEN  -1 

3870  MOVE  0,0 

3880  FOR  1=1  TO  1024  STEP  H 

3890  DRAW  Tl <I),Aw2(I) 

3900  NEXT  I 

3910  PRINT  CHR$<12) 

3920  PRINT  USING  '//////" 

3930  PRINT  ' 

3940  PRINT  ' 

3950  PRINT  ■ 

3960  ON  KEY  4  LABEL  "   RE-PLOT*  GOTO  Replot 

3970  ON  KNOB  .05  GOTO  Vary 

3980  GOTO  Wait 

3990  VaryiPRINT  CHR$(12) 

4000  PEN  -1 

4010  MOVE  0,8 

4020  FOR  1=1  TO  1024  STEP  H 

4038  DRAW  Tl(I),Av2(I) 

4048  NEXT  I 

4058  PEN  1 

4060  C2=KN0BX 

4070  Centerhne=Centerline+C2/120 

4080  MAT  Av2=  Avl/(Centerhne) 

4090  GOTO  3870 


'ERASES  LAST  PLOT 


USE  THE  KNOB  TO  RE-NORMALIZE" 
KEY  *  4  WILL  RE-PLOT  USING" 
UPDATED  SCALE" 


(ERASES  LAST  PLOT 


!KNOB  USED  TO  VARY  NORMALIZING  VALUE 

!NEU  NORMALIZED  DATA  ARRAY 

'PLOTS  THE  NEW  NORMALIZED  ARRAY 


4100  Replot:    ! 

4118    PRINT  CHRSU2) 

4120    CALL  Pl»t(Av2(*>) 

4136    CALL  Result(Av2(*),TK*),H) 

4140    GSTORE  Graf<*) 

4150    GOTO  Chanqe 

4160  Subent: 


i 


4170  SUBEND 

4180  SUB  Menuel 

4198  PRINT  USING 

4200  PRINT 

4210  PRINT 

4220  PRINT 

4230  PRINT 

4240  PRINT 

4250  PRINT 

4268  PRINT 

4270  PRINT 


YOU  CAN  RE-NORMALIZE* 

**************** ********************************* **********  X" 

YOU  CAN  RE-AVERAGE  (New  Data  Only)" 

***************************************** **HX**i**********X" 

PRESS  KEY  *  4  TO  LOOK  AT  OTHER  DATA' 
OLD  OR  HEU  ' 

MOTOR  /  EXHAUST" 
***************************************** ****x****x********x' 
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4?30 
4290 
4300 
4310 
4320 
4330 
4340 
4350 
4360 
4370 
4380 
4390 
4400 
4410 
4420 


*  5 


PRINT  '  PRESS  KEY 

PRINT  '  FOR  THE  '   TWO-ANGLE   '  METHOD" 

PRINT  '  OF  PARTICLE  SIZING'       * 

PRINT  "a******************************** ******** **h********h****«' 

PRINT  '  TO  STORE  THE  REDUCED  DATA  PRESS  KEY  *  6 

PRINT  *************************************************************' 

PRINT  '  TO  END  THIS  SESSION  PRESS  KEY  *  7" 

SUBEND 

SUB  Displayl(Old) 

TO  REDUCE  NEW  DATA  PRESS  KEY  *  1" 


PRINT  ' 
PRINT  " 
PRINT  ' 
PRINT  "   ■ 
ON  KEY  1  LABEL 
ON  KEY  2  LABEL 

4430  Wait :  GOTO  Wait 

4440  New:  PRINT  CHR$( 

4450  PRINT 

4460  PRINT 

4470  PRINT 

4480  PRINT 

4490  PRINT 

4500  PRINT 

4510  PRINT 

4520  PRINT 

4538  Old=0 

4545  SUBEXIT 


TO  REVIEW  PREVIOUSLY  REDUCED  DATA  PRESS  KEY  *  2' 

NEW  DATA"  GOTO  New 
OLD  DATA"  GOTO  Renew 


12) 


PUT  THE  DISK  IN  THE  LEFT  DRIVE  IF  IT  IS  NOT  ALREADY' 


ENTER  THE  NAMES  OF  THE  FIRST  AND  SECOND  FILES,' 


EACH  FILE  HAS  DATA  FROM  BOTH  DIODE  ARRAYS. 


FIRST  IS  NO-PARTICLES  &  NEXT  IS  PARTICLES DH,D2$" 


4556  Review, 
4560 


01d=l 


4578 
4580 
4590 
4600 
4610 
4626 
4630 
4640 
4650 
4668 
4470 
4680 
4690 
4700 
4710 
4728 
4730 
4740 
4750 
4768 
4770 


'THIS  VARIABLE  IS  PASSED  TO  THE  MAIN  PROGRAM  TO  INDICATE 
!THAT  THE  DATA  TO  BE  READ  IS  ONE  (1)  SCAN  AND  THAT 
!THE  AVERAGING  ROUTINES  UILL  NOT  BE  USED 
SUBEND 

SUB  Readata 

COH  /Readata/  B,P,H,Q3$[20],Q4$[20],Zz$[20],YK*)  BUFFER,Y2(»)  BUFFER 


Xyz=l 

B=8 

R  1=65536 

IF  P=2  THEN  Xyz=4097 

IF  P=2  THEN  R 1=32768 

IF  P=2  THEN  B=4 

PRINT  "  ■ 

PRINT  ■ 

ASSIGN  gDiskl  TO  Q3$iZz$ 

ASSIGN  ?Disk2  TO  94$iZz$ 

ASSIGN  BBuffl  TO  BUFFER  Yl(«) 

ASSIGN  gBuff2  TO  BUFFER  Y2(«) 

CONTROL  ?Diskl,5;Xyz 

TRANSFER  SDiskl  TO  §Buffl;COUNT  Rl 

WAIT  FOR  EOT  gDiskl 

CONTROL  gDisk2,5;Xyz 


!  FILE  POINTER  VARIABLE 

!  NUMBER  OF  SCANS  IN  MOTOR  DATA 

•NUMBER  OF  BYTES      OF  MOTOR  DATA 

! RECORD  *  WHERE  EXHAUST  DATA  BEGINS 

! NUMBER  OF  BYTES   OF  EXHAUST  DATA 

(NUMBER  OF  SCANS  IN  EXHAUST  DATA 


READING  DATA  FROM  FILE  ON  DISK" 


(OPEN  I/O  PATHS 


(READS  NO-PARTICLE  DATA 
!  NOT  AN  OVERLAPPING  TRANSFER 
!Xyz  SETS  DISK  FILE  POINTER  TO 
(EITHER  MOTOR  OR  EXHAUST  DATA 
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4783     TRANSFER  §Disk2  TO  ?Buff2;CCUNT  Rl  IREADS  PARTICLE  DATA 

4799     WAIT  FOR  EOT  *Disk2 

4390     ASSIGN  SDiskl  TO  * 

4810     ASSIGN  gDisk2  TO  *  "JUST  GOOD  PRACTICE  TO  CLOSE 

4820     ASSIGN  BBuffl  TO  *  !I/0  PATHS 

4830     ASSIGN  ?Buff2  TO  * 

4840  SUBEND 

4850  SUB  Plot2(E,D(*),X) 

4860  DIM  Tillel$[25],ThetaSI20] 

4870  CINIT 

4880  IF  X=l  THEN  PLOTTER  IS  705,'HPGL' 

4890  DEG  IDEGREES  fir   LABEL  DIRECTION 

4900  GRAPHICS  ON 

4910  VIEWPORT  35,125,30,85 

4920  Max=10*INT((MAX(D(*))+10)/10) 

4930  WINDOW  D(l, 2), D(E,2),0,Max 

4940  G=(D(E,2)-D(1,2))/(E-1>*4        !AN  X  GRID  LINE  EVERY  4th  POINT 

4950  F=(INT(E/4)-l)*4 

4960  IF  F=0  THEN  F=4 

4970  CLIP  D(4,2)-G,D(F,2)+2*G,0,Max     IMAXES  GRID  UNIFORM 

4980  GRID  G,5,»(4,2)-G,8 

4990  CLIP  OFF 

5000  LORG  8 

5010  LDIR  90 

5028  CSIZE  4, .5 

5030  FOR  1=4  TO  E  STEP  4  IPUTS  NUMBERS  ON  X  AXIS 

5040        MOVE  D(I,2),0 

5050       LABEL  USING  '.DDDD";D(I,2) 

5060  NEXT  I 

5070  LDIR  0 

5080  LORG  8 

5090  FOR  1=10  TO  Max  STEP  10  ! NUMBER  Y  AXIS 

5100       MOVE  D(4,2)-G,I 

5118       LABEL  USING  »4,K"jI 

5120  NEXT  I 

5130  CSIZE  6,. 6 

5140  TitleJ='TWO-ANGLE  METHOD" 

5150  Titlel$="For  Various  Angle  Ratios"  'STRINGS  FOR  PLOTS 

5160  Sizet="SIZE  («icrons)* 

5170  Thta$=*THETA  (rad)" 

5180  Sub$=T 

5190  LDIR  90 

5200  LORG  5 

5210  B=D(4,2)-G-(D(E,2)-D<1,2))/18 

5220  MOVE  B,Max/2 

5230  LABEL  Size$ 

5240  LDIR  6 

5258  A=(D(E,2)+DU,2))/2 

5260  MOVE  A, -Max/4 

5270  LABEL  Thta* 
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5230    LORG  3 

5290    MOVE  A,-tei/4 

5380    LABEL  Sub$ 

5310    LORG  5 

5320    MOVE  A,Ha:*1.15 

5330    LABEL  Title* 

5340    CSIZE  4.5 

5350    MOVE  A,Max*1.05 

53*0    LABa  Titlelt 

5370    PENUP 

5380  SUBEND 

5390  SUB  Distribution(D(»),Tratio,E) 

5400    LINE  TYPE  1 

5410    CLIP  ON 

5420    Tratio*=VAL$(Tratio)  !  ANGLE  RATIO 

5430    HOVE  D(1,2),D(1,1) 

5440    FOR  1=1  TO  E  !THIS  SUBROUTINE  PLOTS  THE 

5450       DRAW  D(I,2),D(I,1)  (PARTICLE  SIZES  DERIVED  USING 

5460    NEXT  I  (VARIOUS  ANGLE  RATIOS  APPLIED 

5470    LORG  4  !TO  THE  DATA  OVER  A  RANGE  OF 

5480    LINE  TYPE  1  (ANGLES 

5490    CSIZE  4, .3 

5530    MOVE  D(H,2),D(I-l,l)+2 

5510    LABEL  TratioS 

5520    PENUP 

5530  SUBEND 

5540  SUB  Tuoangle 

5550  !***************************************************************** 

5560  !***************    SUBPROGRAM  'TUOANGLE'    **********x******** 

5570  !***************     PARTICLE  SIZING  BY     ******************* 

5580  !***************     MEASURING  THE  RATIO     a****************** 

5590  !***************      OF  INTENSITY  AT      ******************* 

5600  !***************       TWO  ANGLES        ******************* 

5610  !***************************************************************** 

5620  OPTION  BASE  1 

5630  COM  /Two/  Avl(1024),M,Ml,F 

5640  COH  /Gauss/  T1(1024),G(1024),L  !THE  GAUSS  IS  NOT  USED  HERE 

5650  !BUT  THE  COM  BLOCK  HAS  THETA 

5660  DIM  D(200,2)  !AND  WAVELENGTH 

5670  PRINT  CHR$(12) 

5680  PRINT  USING  '///////' 

5690  PRINT  '    THIS  SUBPROGRAM  USES  THE  TWO-ANGLE  METHOD  DESCRIBED  BY  BUCHELE" 

5708  PRINT  '  " 

5710  PRINT  '    TO  CALCULATE  PATICLE  SIZE  FOR  VARIOUS  ANGLE  RATIOS  AND  ANGLES.1 

5720  PRINT  '  ' 

5730  PRINT  '        IT  IS  HOPED  THAT  THE  CURVES  WHICH  RESULT  WILL  SHOW  ' 

5740  PRINT  '  ■ 

5750  PRINT  "  WHICH  SIZE  IS  THE  MOST  PROBABLE.' 

5760  PRINT  *  ' 

5770  PRINT  "AFTER  NOTING  FROM  THE  RAW  DATA  WHICH  ANGLES  CONTAIN  THE  CENTER  LOBE" 
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5780  PRINT  " 

57<?Q  PRINT  '  ENTER*************THE  SMALLEST  USEABLE  ANCLE," 

S800  PRINT  ■  THE  SMALLEST  ANGLE  RATIO, AND" 

5810  PRINT  '  THE  STEP  SIZE  BETWEEN  ANGLE  RAIOS' 

5820  PRINT  "  YOU  WISH  TO  EXPLORE" 

5830  PRINT  "         EXAMPLE*************  ,012, 1.2, .4' 

5840  INPUT  '  ENTER         thetal,theta-ratio,theta-step",Q,A,B 

5850    OFF  KEY 

5860    X=0  !GRAPH  WILL  APPEAR  ON  SCREEN 

5870  Begin:  "RUNNING  CONTINUES  HERE  WO  A  HARD  COPY  IS  DESIRED 

5880       (Van  De  Hulst  and  Guwprech  &  Sleepevich  explain  that  the  chanae  in 

5890       'wavelength  of  the  bean  is  accounted  for  by  dividing  b?  the  index 

5900        !of  refraction  of  the  nediun,  *********L=L/M1  ********** 

5910  C=(L/Ml/.57/PI)*2       ! see  page  15  of  nasa  tech  paper  2156 

5920  FOR  N=l  TO  1024  Its  see  this  is  a  convenient  constant 

5930     IF  Q(TKN)  THEN  5950  (FINDS  POSITION  OF  KINIHUN  ANIGE 

5940  NEXT  N 

5950  FOR  Tratio=A  TO  3  STEP  B  (VARIOUS  ANGLE  RATIOS 

5960    FOR  J=N  TO  1010/Tratio  STEP  10  !SETS  THE  RANGE  OF  POSSIBLE 

5970  !  ANGLES  (THETA  1) 

5980       Thtal=TKJ) 

5990       Thta2=TKJ*Tratio)  !  THETA  2 

6000       Deltatheta=Thta2A2-ThtalA2 

6010       11=0 

6020       12=0 

6030  FOR  I=J-5  TO  J+5  ITHESE  2  LOOPS  AVERAGE  A  FEW 

6040  Il=Il+Avl(I)  'INTENSITIES  IN  THE  HOPE 

6050  NEXT  I  !OF  A  STEADIER  CALCULATION 

6060  11=11/11 

6070  FOR  I=INT(J*Tratio-5)  TO  INT(J*Tratio+5) 

6080  I2=I2+AvlfI>  !  DONE  HERE  FOR  THETA2  FOR 

6090  NEXT  I  !  THE  GIVEN  ANGLE  RATIO 

6100  12=12/11 

6110  E=(J-N)/10+1  iTHIS  IS  A  COUNTER  FOR  THE  ARRAY  CONTAINING 

6120  (PARTICLE  SIZE  AND  THETA1  FOR  THE  GIVEN 

6138  (ANGLE  RATIO 

6140  IF  11(12  OR  11=0  OR  I2<=0  THEN  D(E,1)=0   (ALLOWANCE  FOR  IF  THE  DATA 

6150  IF  11(12  OR  11=0  OR  I2(=0  THEN  GOTO  Xcowp     IIS  NOT  WELL  BEHAVED 

6160  D(E,l)=SGR(-C/Deltatheta*L0G(I2/Il))  ! TWO  ANGLE  METHOD 

6170  (VALUE  OF  DIAMETER  based  an  INTENSITY  RATIO 

6188  (For  a  given  ANGLE  RATIO  and  ANGLE  THETA1 

6190  Xconp:      D(E,2)=Thtal 

6200       Spar=PI*D(E,l)*Ml/L    (THIS  IS  PI*D/LA"MBDA  -  THE  SIZE  PARAMETER 

6210       Tbar=Spar*Thta2  (THETA  BAR  FOR  THE  LARGER  ANGLE 

6220       IF  Tbar>3  THEN  J=1010/Tratio      (THIS  ENDS  THE  DO  LOOP  FOR  THIS 

6230       IF  E=l  AND  Tbar)3  THEN  6330       (ANGLE  RATIO  SINCE  THE  GAUSSIAN 

6240  (the  above  line  ends  all  calculation  if     (IS  NOT  VALID  WHEN  Tbar  >  3 

6250  (the  first  elenent  (E=l)failed  the  test 

6260    NEXT  J 

6270  !   FIRST  TIME  THROUGH-  D  HAS  THE  MOST  ELEMENTS  IT  WILL  HAVE 
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6290  !   £0  SET  UP  THE  GRAPH  USING  D's  PRESENT  PARAMETERS 

6290  PRINT  CHR$(12) 

6383    IF  Tratio=A  THEN  CALL  Plot2(E,D(») ,X) 

6310    CALL  Distribution<D<*),Trau»,E) 

6320  NEXT  Tratio 

6330    CN  KEY  2  LABEL  'HARD  COPY"  GOTO  Hard 

6340    ON  KEY  3  LABEL  "HENUE  '  GOTO  Sebexit 

6350  PRINT  USING  '/////' 

6360  PRINT  "  YOU  CAN  GET  A  HARD  COPY  BY  PRESSING  KEY  *  2' 

6370  PRINT  '  OR  EXIT  THIS  ROUTINE  BY  PRESSING  KEY  t  3" 

6380  Standby -.GOTO  Standby 

6390  Hard:X=l  !A  VARIABLE  TO  CONTROL  PLOTTER 

6400     GOTO  Begin 

6410  Subexit:  GINIT 

6420       GRAPHICS  OFF 

6430  SUBEND 

6440  SUB  Shift 

6450  COM  /Readata/  B,P,H,Q3$[20],94$[20],Zi$[20],YH*)  BUFFER, Y2(»)  BUFFER 

6460  DIM  E(l:9192) 

6470  PRINT  CHRS(12) 

6480  PRINT  USING  "///////' 

6490  PRINT  '  RAW  DATA  IS  BEING  SHIFTED  TO  CORRECT  FOR  SMALL  GAPS  BETWEEN  SCANS" 

6500  PRINT  '" 

6510  PRINT  ■  BE  PATIENT,  IT'S  A  LONG  SET  OF  DO-LOOPS" 

6520      !      THERE  ARE  SOME  SHALL  GAPS  BETWEEN  SCANS  AND  THIS  SUBROUTINE 

6530      !      SHIFTS  THE  DATA  SO  THAT  THE  FIRST  DIODE  DATA  POINT  IS  MOVED 

6540      !      TO  THE  VERY  BEGINNING  OF  ITS  1624  BLOCK  IN  THE  OVERALL  ARRAY 

6550      !      .THE  FIRST  SET  IS  RIGHT  ON,  THE  NEXT  IS  ONE  OFF, THE  THIRD  IS 

6560      !      TWO  OFF  ,,SO  FORTH.  THIS  MAY  NOT  MATTER  WITH  OUR  RESOLUTION 

6570      !      AND  IS  PROBABLY  DUE  TO  THE  MEMORY  CARD  CYCLING  AT  THE  END  OF 

6580      !      SOME  SCANS.  SEE  THE  CIRCUIT  TIMING  DIAGRAM  IN  THE  THESIS. 

6590     SELECT  B 

6600     CASE  4  ! EXHAUST  DATA  HAS  4  SCANS 

6610         H=3 

6620     CASE  8  (MOTOR  DATA  HAS  8  SCANS 

6630         M=7 

6640     END  SELECT 

6650   FOR  K=l  TO  2 

6660   IF  K=l  THEN  MAT  E=  Yl      IOPPERATES  ON  NO-PARTICLE  AND  PARTICLE  SETS 

6670   IF  K=2  THEN  HAT  E=  Y2 

6680         FOR  J=0  TO  M 

6690         IF  M=7  AND  J<=3  THEN  6800  !One  4096  Block  Doesn't  Need  Shifting 

6700  FOR  I=(J)»1024+1  TO  (J+l)*1024         ! BLOCKS  OF  1024 

6710  IF  M=3  THEN  L=I+J     'Array  B  Has  the  worst 

6720  IProblert  With  Shifting  Data 

6730  IF  M=7  THEN  L=I+1      lArray  D  is  Allways  off  by  one 

674J  IF  LXJ+1X1024  THEN  L=(J+1)*1024 

6750  ! JUST  TO  AVOID  PROGRAMMING  ERROR  AT  THE  END  OF  THE  ARRAY 

6760  E(I)=E(L)  ITHIS  SHIFTS  THE  DATA 

6770  "DEPENDING  ON  'L', ARRIVED  AT 
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6780  !BY  OBSERVING  RAW  DATA 

6790  NEXT  I 

6800         NEXT  J 

6810     IF  K=l  THEN  MAT  Yl=  E 

6820     IF  K=2  THEN  HAT  Y2=  E 

6830  NEXT  K 

6840  SUBEND 


6850  SUB  Store 

6860 

6870 

6880 

6890 

6900 

6910 

6920 


THE  LARGE  ARRAY  CONTAINING  MULTIPLE  SCANS  HAS  BEEN  REDUCED 
TO  A  MEAN  SCATTERING  PROFILE  BY  AVERAGING  AND  FILTERING.  IF 
YOU  FEEL  THAT  STORING  THIS  DATA  IS  NECESSARY  THIS  ROUTINE 
DOES  IT,  IT  SAVES  DISK  SPACE  TO  STORE  THE  RESULTS  THEN 
ELIMINATE  THE  RAW  DATA  IF  YOU  FEEL  CONFIDENT  THAT  THE 
REDUCTION  IS  THE  BEST  THAT  CAN  BE.  IN  OTHER  WORDS, 
DO  NOT  PURGE  A  RAU  DATA  FILE  UNLESS  YOU  ARE  ABSOLUTELY  SURE 
YOU  WON'T  WANT  TO  REDUCE  IT  AGAIN. 

6930  COM  /Hrdqaus/  Av2(»)  !Av2<»)  IS  THE  REDUCED  DATA 

6940  DIM  DataU : 1024)  BUFFER  'BUFFER  USED  TO  TRANSFER  TO  DISK 

6950  MAT  Data=  Av2 

6960  PRINT  CHR$(12) 

6970  PRINT  USING  *//////" 

6980  PRINT  *   A  SUGGESTED  HETHOD  OF  NAMING  REDUCED  DATA  FILES  IS  AS  FOLLOWS" 

6990  PRINT  '  M MOTOR  BEAM" 

7000  PRINT  '  X EXHAUST  BEAM" 

7010  PRINT  '  C CALIBRATION  ,  IF  NO  'C  THEN  AN  ACTUAL  FIRING  IS  ASSUMED' 

7020  PRINT  '  G GLASS  BEADS  ' 

7030  PRINT  *  IF  NO  'C  THEN  'G'  STANDS  FOR  'GAP'  PROPELLANT" 

7043  PRINT  '  A ALUMINUM  OXIDE" 

7050  PRINT  '  PPP CHAMBER  PRESSURE  FOR  RUN  OR  OTHER. .. .(NOZZLE  TYPE)" 

7060  PRINT  '       PARTICLE  SIZE  FOR  CALIBRATION" 

7070  PRINT  "  MHDD- — MONTH, DAY  OF  RUN  OR  CALIBRATION" 

7080  PRINT  " 

7090  PRINT  '  EXAMPLE:   MCA125JN12" 

7100  PRINT  "" 

7118  PRINT  '  MOTOR  BEAM  CALIBRATION  USING  ALUMINUM  OXIDE  FROM  1  TO  25  MICRONS" 

7120  PRINT  "  ON  JUNE  12" 

7130.  PRINT  ■  EXAMPLE:   XG550AU10" 

7140  PRINT  " 

7150  PRINT  "  EXHAUST  BEAM  DATA  FOR  GAP  PROPELLANT  AT  550  psi  ON  AUGUST  10" 

7160  PRINT  *  PLACE  A  DISK  IN  THE  RIGHT  HAND  DRIVE' 

7170  PRINT  '  ENTER  THE  FILENAME  YOU  WISH  TO  USE  FOR  THIS  REDUCED  DATA" 

7180  INPUT  A* 

7198  CREATE  BDAT  AS, 512, 16 

7200  ASSIGN  BDisk  TO  AS 

7210  ASSIGN  BBuff  TO  BUFFER  DataU) 

7220  CONTROL  BBuff ,3;1, 8192,1  IBUFFER  IS  FULL 

7230  TRANSFER  BBuff  TO  BDisk ;COUNT  8192 

7240  WAIT  FOR  EOT  BDisk 

7250  ASSIGN  BBuff  TO  * 

7260  ASSIGN  BDisk  TO  * 

7270  SUBEND 


83 


7280 
7290 
7300 
7310 
7320 
7330 
7340 
7350 
7360 
7378 
7380 
7390 
7400 
7410 
7420 
7430 
7440 
7450 
7460 
7470 
7480 
7490 
7588 
7518 
7528 
7530 
7540 
7550 
7560 
7578 
7580 
7598 
7600 
7618 
7620 
7638 
7640 
7658 
7660 
7670 
7680 
7690 
7780 


THE  DISK  WITH  THE  REDUCED  FILE  SHOULD  BE  IN  THE  RIGHT-HAND  DRIVE" 

ENTER  THE  FILENAME  OF  THE  REDUCED  DATA' 
TO  BE  REVIEWED' 


!THIS  IS  AN  EMPTY  BUFFER 
!1024*3=NUMBER  OF  BYTES 


SUB  Revieu(Avl(*)) 

DIM  Data(l:1024)  BUFFER 

PRINT  '  EACH  REDUCED  FILE  CONTAINS  ONE  SCAN  ONLY.  YOU  MUST  KNOU  IF  IT1 

PRINT  '     IS  EXHAUST  ,  MOTOR  CAVITY  ,  OR  CALIBRATION  DATA." 

PRINT  " 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT  ■* 

INPUT  A$ 

ASSIGN  gDisk  TO  At 

ASSIGN  BBuff  TO  BUFFER  Data!*) 

CONTROL  8Disk,5;l 

CONTROL  gBuff ,3;1 ,0 ,1 

TRANSFER  SDisk  TO  PBuff jCOUNT  8192 

WAIT  FDR  EOT  gDisk 

ASSIGN  gDisk  TO  * 

ASSIGN  gBuff  TO  * 

MAT  Avl=  Data 

SUBEND 

SUB  Filter(A(*),Fil) 

DIM  E(l:1024) 

PRINT  CHR$(12) 

PRINT  USING  '////////' 

PRINT  • 

FOR  J=1  TO  Fil 

FOR  1=1  TO  1014 

B=I+1 

C=I+2 

D=I+3 

F=I+4 

G=I+5 

H=I+6 

K=I+7 

L=I+8 

N=I+9 

P=I+10 


FILTERING  THE  SCATTERING  PROFILE" 


!THIS  IS  A  Svnetric  Hovinq  Averaqe 
!TYPE  OF  DIGITAL  FILTER,  EACH 
[DATA  POINT  IS  EQUALLY  WEIGHTED 
!IN  THIS  CASE  BUT  THIS  CAN  BE 
(CHANGED  IF  ONE  DETERMINES  THAT 
•FEWER  POINTS  WITH  UNEQUAL  WEIGHTS 
[WOULD  BE  FASTER  OR  GIVE  BETTER 
[RESULTS,  THIS  TYPE  OF  FILTER  WAS 
[USED  SINCE  IT  INTRODUCES  NO  PHASE 
[LAG      (Anqular  Resolution) , 


E(G)=(A(I)+A(B)+A(C)+A(D)+A(F)+A(G)+A(H)+A(K)+A(L)+A(N)+A(P))/11 
NEXT  I 
HAT  A=  E 
NEXT  J 
SUBEND 
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